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STUDY  OP  THE  OPERATION  OF  A  CENTRIFUGAL  COMPRESSOR  WHEEL* 


S.  A.  Anisimov 

The  problem  of  creating  highly  efficient  centrifugal  compressors, 
at  the  modern  level  of  the  development  of  the  theory  of  these  machines. 
When  theoretical  methods  alone  do  not  assure  sufficiently  reliable 
results,  can  be  successfully  solved  only  by  means  of  appropriate 
experimental  work.  Of  great  significance  along  these  lines  are 
experimental  studies  of  centrifugal  compressor  wheels;  when  evaluating 
the  efficiency  of  operation  of  the  wheels  it  is  necessary  to  take  into 
account  not  only  the  basic  characteristics  of  the  wheel  itself,  but 
also  the  effect  that  the  wheel  has  on  the  operation  of  the  fixed 
elements,  mainly  on  the  diffusor.  In  other  words,  the  operational 
efficiency  of  a  wheel  should  be  Judged  not  only  on  the  basis  of  its 
characteristics,  but  also  on  the  basis  of  the  structure  of  the  flow 
behind  the  wheel. 

The  studies  were  made  on  the  ETsK-1  stand  at  the  Compressor 
Laboratory  of  the  Kalinin  Leningrad  Polytechnic  Institute  [l].  In 
all,  we  studied  10  double-shrouded  wheels  with  single-stage  cascades; 

i 

*  Docent  K.  P.  Seleznev  and  Senior  Engineer  F.  S.  Rekstin 
participated  in  the  work. 
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the  only  difference  was  in  the  number  of  blades.  Figure  1  gives  the 
basic  geometric  parameters  of  these  wheels.  In  the  article  we  give 
experimental  data  on  a  single  wheel  (z  =  =  zg  «=  16),*  since  the 

principles  established  for  the  other  wheels  are  principally  of  a 
similar  nature. 

Special  observations  were  made  for  cases  when  the  nature  of  the 
curves  changed  somewhat  with  a  change  in  the  number  of  blades  z.  The 
question  of  the  influence  of  the  number  of  blades  on  the  efficiency 
of  the  wheels  is  treated  separately  in  detail  in  the  second  article 
in  this  series  (see  p.  20). 

A  diagram  of  the  blading  of  the  ETsK-1  test  stand  and  the  control 
sections  for  measuring  the  individual  flow  parameters  are  given  in 
Fig.  2  and  have  been  treated  in  detail  by  Rekstin  [i]. 

At  cross  section  A-A  (atmospheric  conditions)  we  measured  the 
temperature  T  with  precision  thermometers  graduated  to  0.1°  C,  and 
the  barometric  pressure  P  .  At  cross  section  0-0,  at  the  wheel  eye, 
we  measured  the  total  pressure  with  a  cylindrical  ram-air  pipe,  and 
the  static  pressure  through  holes  0.7  mm  in  diameter  in  the  housing 
of  the  inlet  duct  and  the  fairing.  At  cross  section  3-5,  behind  the 
impeller,  we  determined  the  fields  of  temperature,  angles,  and  total 
and  static  pressure  using  temperature  probes  and  pneumometric  heads. 

In  addition  we  measured  the  static  pressure  through  holes  0.7  mm  in 
diameter  in  the  forward  wall  of  the  diffuser  and  in  the  diaphragm. 

The  air  output  was  determined  using  interchangeable  metering 
nozzles,  prepared  and  installed  in  the  induction  manifold  according 
to  the  existing  norms.  The  internal  power  was  determined  by  two 

*  This  wheel,  with  z  =  16,  was  the  best,  from  the  standpoint  of 
the  realizable  efficiency  and  head,  of  all  the  wheels  investigated. 
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methods:  by  the  drop  in  total  pressure  and  by  the  reactive  moment  in 
the  housing  of  the  electric  .motor  with  a  pivoting  stator. 


The  number  of  revolutions  of  the  shaft  of  the  experimental  stage 
was  measured  from  the  readings  of  an  electrical  tachometer  and  a 
total  revolution  counter. 

The  wheels  were  tested  for  a  constant  number  of  revolutions, 

n  =  14,900  rpm  (u2  =  214*5  m/sec). 

•  To  evaluate  the  efficiency  of  the  wheels  we  used  the  following 

criterional  parameters:  the  adiabatic  efficiency,  by  the  total  para- 

— * 

meters  ti*v  and  its  relative  value  vj  the  head  coefficients 
and  <p;  the  degrees  of  reaction  Cl  and  p ;  the  coefficients  of  relative 
nonuniformity  of  flow  kQ,  k^,  k£  and  k^  ;  the  absolute  nonuniformity 
of  flow  Aa;  and  the  M  and  Re  numbers.  As  the  independent  determining 
parameter  we  selected  the  flow  factor  q>p  (sometimes  cpp)  and  the 
angle  of  attack  i^  which  is  practically  always  associated  with  it.* 

In  our  studies  the  Re  number  was  arbitrarily  determined,  as  is 
often  done  in  the  theory  of  centrifugal  compressors,  according  to  the 
formula 

Rp_“«d.  _  .  . 


and  was  within  the  limits  (3.I-5.6)  x  10^. 

The  parameters  were  averaged  in  various  ways.  The  average  values 
of  the  parameters  PQ,  P,  and  Tq  were  determined  as  the  arithmetic 
means.  The  average  values  of  the  parameters  f  and  h&^  were  found 


*  The  arbitrary  designations  used  in  the  article,  with  few  excep¬ 
tions,  are  the  same  as  the  generally-accepted  ones  [2].  If  not, 
explanations  are  given  in  the  text. 

The  subscript  of  the  parameters  Indicates  the  number  of  the 
control  section  to  which  the  parameters  pertain.  The  parameters 
behind  the  wheel,  in  cross  section  5-5,  are  written  without  a  sub¬ 
script. 
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from  previously  averaged  values  of  PQ,  P,  and  Tq.  The  average  values 
of  the  parameters  t}*^,  f*,  9,  a,  h*^,  h^,  and  0.  were  found  by  averag¬ 
ing  their  local  values  with  respect  to  the  flow.  '  The  average, flow 
values  Q  and  9p  were  determined  by  measuring  with  a  planimeter. 

1.  Results  or  the  Investigation 


1. 

where 


Figure  3  shows  the  dependences  t)*^  =  f(<pp  )  and  r)a^  =  f(<pp  ) 


Tp<> 


Nad  Cihad  ha» 

G'h,  hi  ' 

(2) 

—  •  » 

r  in!  n 

(2») 

(3) 

(4) 

(5) 


For  the  value  of  9  when  T)*-.  and  we  used,  respectively, 

0.0  ITlaX  luciX 

the  designations  9pQOpt  and  9^. 


-*  .  . 

The  curve  ^  ~  f" C 'Pp  )  bas  an  extreme  value.  The  drop  in  effi¬ 
ciency  to  the.  left  and  to  the  right  of  the  extreme  value  is  explained 
by  the  large  positive  and  negative  angles  of  attack  which  cause  a 
deterioration  in  the  aerodynamics  of  the  flow  in  the  wheel,  and  by 
the  appearance  of  developed  separation  zones. 


-* 


The  experimental  data  on  t]  ^  obtained  by  various  methods  —  accord 
ing  to  temperatures  (t^  t )  and  using  a  pendulum-motor  (t^  )  —  are 

in  good  agreement.  The  maximum  value  of  the  efficiency  attained  for 


the  given  wheel  t]*^  t  _  =  0.  9. 


max 


The  nature  of  the  curve  T]*^ 
blades.  With  a  decrease  in  z  the  maximum  efficiency  is  weakly 


f(9  )  depends  on  the  number  of 

p0 
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expressed,  while  when  z  2  6  it  is  generally  not  expressed  at  all. 

2.  The  head  coefficients  shown  in  Fig.  3  were  determined  from 
the  following  formulas: 


where 


-H 


(6) 

(7) 

(8) 

(9) 


The  curves  of  the  head  coefficients  if/*,  ip,  and  cp  have  an  extreme 
value . 

In  the  optimum  regime  (at  T]*^  max)  these  coefficients  were  ip*  = 

=  1.28,  ip  =  0.82,  and  9  =  O.695.  For  the  regime  i1  =  0  we  obtained 

ip*  =  1.15,  ip  =  O.75,  and  9  =  O.615. 

The  decrease  in  the  head  coefficients  ip*  and  9  with  Increased 

flow  in  the  region  9  >  9  can  be  explained  from  an  examination 

Pq  Pocr 

of  the  Stodola  formula  for  the  coefficient  92  (considering  that  the 

ratio  9  /9  is,  for  all  intents  and  purposes,  approximately  constant) 
P2  Pq 


<P2  =  —  =  1 
M- 


n 


sin  P*i  C  9p»  cot  —  I 


n  .  a 
•— smp,,- 


*/>, 


cot  P.>,  •  <Pp,  (10) 


and  on  the  basis  of  the  connection  between  the  coefficients  9  and  ip*: 

r  =  2cp  <l  +  +  pmp)  c.  (11) 


Formula  (ll)  was  obtained  after  transformations  from  Formulas 
(2)  and  (6),  taking  into  account  the  relations 

hi  =  hT  (1  +  +  M- 
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(12) 

(1?) 


where,  strictly  speaking,  by  f*  and  cp  in  Formula  (ll)  we  should  under¬ 
stand  ijj*  and  <p2> 

In  addition,  a  drop  in  and  cp  for  large  ©  (i,  <  0)  is  condi- 

Pq  x 

tioned  by  a  deterioration  in  the  flow  structure  (the  development  of 


separation  zones  in  the  wheel  channels). 

The  decrease  in  the  head  coefficients  for  low  flow-rates  (cp_  < 

pO 

<  <p„  is  explained  by  the  sharp  deterioration  of  the  flow  structure 


P0cr 

in  the  wheel  with  large  positive  angles  of  attack.  In  this  region 
Formula  (10)  cannot  be  used. 

The  curves  ij^  =  f(<p  )  and  q>  =  f ( ep  )  have  approximately  coin- 

p0  p0 

cident  maxima,  while  the  curve  f*  =  f(<p  )  has  great  curvature.  This 

p0 

can  be  explained  by  the  structure  of  Formula  (ll),  since  the  term 
(l  +  Pnp  +  whlch  characterizes  the  hydraulic  energy  losses 

is  always  less  than  unity,  and  in  the  investigated  regime  range  is 
practically  greater  than  0. 5. 

The  curvature  of  the  curves  cp  =  f(q>  )  and  rj/*  =  f(<P  )  depends 

PO  p0 

on  the  number  of  blades.  With  decreasing  z  the  curves  become  more 
sloping,  so  that  at  small  z  the  maximum  of  the  head  coefficients  is 
weakly  expressed,  while  when  z  =  6  it  is  practically  nonexistent. 

3.  The  degrees  of  reaction  fi  and  p  were  determined  from  the 
following  formulas: 


r>  >>T  —  I’d  _  1  h0  . 
M - r-  1  hT> 


Q  - 


had 


hr 
h'ad  —  h<> 


nad 


where 


/'a  — 


e*  —  4 
2  e 


hd  _  »  . 
’ 


=  (<P*+ 


(14) 

(15) 

(16) 


Formulas  (14)  and  (15),  after  appropriate  transformations,  can  be 
reduced  to  the  following  approximate  form  which  is  convenient  for 
analysis : 
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Q  =  1 


yi  +  Tp-Tp, 


2<p 


— 


9*  +  9p  -  9?. 


29  ( 1  +  Pnp  +  P mp)  TIJa 


2(1+  Pnp  +  P  mp)  'H3a 


(14 ») 
(15' ) 


2  2  2 

since,  for  all  practical  purposes,  9  »  9  -  9  ,  and  therefore  the 

Po 

difference  (9  -  9  )  in  these  formulas  can  be  disregarded. 

P  P0 

The  curves  Cl  =  f(cp  )  and  p  =  f(9  )  throughout  the  entire  flow 

Po  p0 

range,  except  for  the  region  of  large  $  ,  are  of  identical  nature. 

p0 

In  this  case,  in  the  region  of  operating  regimes  (when  9  >9  ) 

P0  P0cr 

with  Increasing  9  the  value  of  Cl  continually  Increases,  while  the 
°0 

value  of  p  first  increases  and  then  (at  large  9  )  begins  to  decrease 

?0 

somewhat.  This  latter  fact  can  be  explained  by  the  influence  of 

energy  losses.  Quantitatively,  Cl  and  p  differ  little  from  one  another, 

except  for  the  region  of  large  9  . 

p0 

The  nature  of  the  change  In  the  degree  of  reaction  with  a  change 

In  flow  rate  can  be  explained  using  Formulas  (141)  and  (15').  Since 

(9„  -  9_  )  «  9  ,  and  it  changes  relatively  insignificantly,  we  can 
p  p0 

conclude  that  the  nature  of  the  change  in  Cl  will  be  approximately 
the  reverse  of  the  change  In  the  value  of  9.  The  value  of  p,  a3  can 
be  seen  from  Formula  (15’),  also  depends  on  the  energy  losses;  there¬ 
fore  its  relation  with  9  will  not  be  quite  so  simple.  Since  the 
value  of  (l  +  Pnp  +  Pmp)  Ti*^  is  always  less  than  unity  in  all  operat¬ 
ing  regimes,  and  since  the  adiabatic  efficiency  takes  into  account, 
in  addition  to  leakage  and  web  losses,  losses  in  the  blading  of  the 
wheel,  it  follows  from  Formulas  (14')  and  (15’)  that  Cl  >  p.  This  is 
also  confirmed  by  experimental  data:  in  the  optimum  regime  Cl  =  0.68 
and  p  «  0.64,  while  In  the  regime  1^  =  0  we  obtained  Cl  =  O..75  and 

p  *  0.68. 

4.  The  curve  9p  =  c^/Ug  =  f(<Pp  )  has  a  dependence  which  is 


close  to  linear,  with  a  slight  downward  convexity  (Fig.  2).  In 
principle,  this  curve  should  pass  through  the  coordinate  origin. 


Quantitatively,  <pp  and  <pp  differ  little  from  one  another,  but  in 

this  case  for  all  operating  regimes  we  have  the  inequality  ®  '.<  <P  . 

P  Po 

Their  relation  in  the  general  case  is  defined  as  follows: 


_  cr  _  .  _£o  _  _Y»^ 

<Pp0  C‘0  C»  F  V  ' 


Analysis  of  the  obtained  experimental  data  shows  that  for  the 

given  wheel  the  ratio  <p_/q>  »  const. 

P  Pq 

5*  The  angle  of  attack  when  the  flow  enters  the  blades  of  the 
wheel  is  determined  from  the  formula 

h  —  fig,  *—  P»  =*  arctmq)^  •  ~ .  ( 17  ) 

In  this  formula  the  inflow  angle  was  determined  without  taking 

into  account  leakage  through  the  packing  on  the  covering  disk  and 

blockage  of  the  channel;  in  addition,  it  was  assumed  that  70  ■»  7^. 

The  curve  i.  =  f  (<p  )  has  slight  downward  convexity.  With 

■L  Pq 

increasing  <p  the  algebraic  value  of  i.  continually  decreases,  which 
p0  1 
also  follows  from  Formula  (17). 

The  maximum  efficiency  was  reached  at  =  +7°.  The  location 
of  the  optimum  regime  in  the  region  of  positive  angles  of  attack  can 
be  explained,  physically,  by  the  fact  that  for  normally  backward  bent 
vanes  having  negative  angles  of  curvature  of  the  radial  cascade  pro- 


file  8  -  Pg2  -  Pg2np,*  positive  1* 


decrease  the  fluid  deflection  in 


*  Pg  np  13  the  geometric  angle  at  which  the  vanes  come  from  the 

wheel,  in  which  case  for  the  given  radial  cascade  we  have  rectilinear 
vanes.  The  values  of  this  angle  for  vanes  along  the  peripheral  arc  is 
determined  from  the  formula  cos  B  =  cos  p  /(D^/D. ).  Since  in 

g2nP  c'  1' 


centrifugal  compressors  Dp/D.  >  1,  then  B  >  B  .  If  D0/D.  =  1, 

g^L  2  1 

we  have,  approximately,  p  =  p  in  the  axial  cascade,  and  conse- 

e>2  P 

quently  the  angle  0  =  p  -  p  . 

S1  gl 


’g2nP 
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the  radial  cascade  e 


P2  -  Pg2np  +  i*  (where  P2  =  arc  tan  <P2p/l  -  <p2). 


which  evidently  has  a  favorable  effect  on  the  aerodynamics  of  the  flow 

=  65°57'  and,  conse- 


in  the  wheel  channels.  For  the  given  wheel  p 
quently,  0  =  49°  -  63°57'  =  -i4°57,« 


g2nP 


The  regime  r)*^  max  has  the  corresponding  angle  P2  =  37  JO*  and, 

consequently,  e  =  37°30»  -  63°57 1  +  7°  «  -19°27*.  Figure  3  shows  the 

dependences  e  =  f(<p  )  and  p„  =  f(cp  ).  Experimental  data  show  that 

Pq  d  Po 

for  the  investigated  series  of  wheels  the  nature  and  position  of  the 

curves  e  =  f(<p  )  and  p_  =  f(qp  )  are  practically  independent  of  the 
PO  d  p0 

number  of  vanes.  For  all  the  investigated  wheels  there  is  a  sharp 
increase  in  the  angle  for  low  flow-rates,  right  up  to  <p  »  0.3,** 

^  Po 

while  the  minimum  values  of  the  absolute  magnitudes  of  the  angles  s 

are  in  the  region  qp  =  0.18-0.23,  close  to  cp  ..  We  should  mention 

P0  PO0?1' 

that  for  the  examined  wheel  with  z  =  16  which,  -as  we  have  already 
stated,  was  optimum  from  the  standpoint  of  efficiency,  the  minimum 
value  of  the  absolute  magnitude  of  e  (|em-j_n!  w  17°)  was  the  smallest 
of  all  the  investigated  wheels. 

6.  The  general  expression  for  determining  the  coefficients  of 

the  relative  nonuniformity  of  flow  according  to  the  angle  k  .  the  flow 

cc 

rate  k~,  and  the  absolute  and  internal  heads  kf  and  k.  ,  shown  in 
G  ^ad  ^i 

Fig.  3,  had  the  form 


kn  —  • 


(18) 


where  II  is  a  measureable  parameter,  averaged  over  the  width  of  the 
channel;  and  IImax  and  are  the  maximum  and  minimum  values  of  the 


*  For  an  axial  cascade  (i.e.,  when  p 


as  6  )  S 

g2nP 


Pi  +  1, 


and:e  -  P2  - 

**  For  other  z's  in  the  region  <p_  >0.3  the  intensity  of  the 

p0 

increase  in  p2  also  decreases,  but  not  quite  so  strongly  as  was  the 
case  with  the  given  wheel. 


-9- 


parameter  n  over  the  width  of  the  channel. 

The  larger  the  coefficient  kjj,  the  greater  the  relative  nonuniform- 

ity. 

The  relative  nonuniformity  of  flow  with  respect  to  energy, 

characterized  by  the  coefficients  k^*^  an<i  ^  respectively  as  a 

function  of  the  regime,  ohanged  but  little.  The  absolute  values  of 

the  relative  nonuniformity  of  the  flow  with  respect  to  energy  are 

very  small  (the  coefficients  k,*  and  k.  are  close  to  zero).  A 

ad  i 

certain  increase  in  k^*  and  k^  is  observed  for  negative  i^  (large 

ad  i 

cp  ),  and  a  very  insignificant  one  is  observed  in  the  region  of  posi- 
P0 

tive  i.  (small  9  ).  Such  an  increase  in  the  coefficients  of  relative 

1  ”0 

nonuniformity  is  evidently  a  result  of  the  deterioration  of  the  flow 
structure  of  the  impeller  for  rather  large  positive  and  negative 
angles  of  attack. 

The  relative  nonuniformity  of  flow  with  respect  to  flow-rate  and 

angle,  characterized  by  the  coefficients  k^  and  ka,  for  angles  of 

attack  i^  «  0  and  regimes  close  to  optimum,  i.e.,  when  t|*^  max  is 

slight  compared  with  other  regimes  (kQ  »  ka  .«  0.25-0.4).  For  large 

positive  angles  of  attack  the  coefficients  of  nonuniformity  k^,  and  kQ 

greatly  increase,  reaching  a  comparatively  high  value  (kQ  ka  «  4-4.5)* 

The  sharp  Increase  in  the  coefficients  kQ  and  kQ  for  low  flow-rates 

is  explained  not  only  by  the  deterioration  in  the  flow  structure,  but 

also  by  the  decrease  in  the  absolute  average  values  of  G  and  a.  For 

negative  angles  of  attack  with  increasing  9  the  relative  nonuniform- 

p0 

ity  with  respect  to  angle  and  flow  rate  increase  insignificantly. 

When  analyzing  the  operation  of  the  stage  with  the  vaned  diffuser, 
a  great  role  is  played  by  the  numerical  values  of  the  deviations  of 
the  angle  from  the  average  angle.  Figure  3  shows  the  dependences  of 


-10- 


a  -  f(®  )  and  the  absolute  nonuniformity  with  respect  to  the  angle 

p0 

Aa  =*  amax  -  ■  f(<Pp  )•  The  graphs  show  that  a,  with  a  change  in 

flow  rate,  changes  within  broad  limits  and  continually  increases  with 

increasing  ©  .  The  curve  Aa  =  f(<p  )  has  a  minimum  for  regimes 

p0  p0 

corresponding  approximately  to  ti*^  max  and  i^  »  0.  The  increase  in 
Aa  in  the  region  of  large  positive  and  negative  angles  of  attack  is 
explained  by  the  deterioration  of  the  flow  structure  within  the 
Impeller  due  to  nonoptimal  flow  past  the  vaned  apparatus  and  the 
influence  of  bending  of  the  flow  at  the  wheel  inlet. 

In  the  region  of  the  investigated  regimes,  the  maximum  value  of 
the  absolute  nonuniformity  of  flow  with  respect  to  the  angle  Aa  reached 
20°,  while  for  the  regimes  max  and  i^  =•  0  it  was  within  limits  of 
70°. 


7.  The  curves  <p  =  f(X)  (Fig.  4)  show  that  for  average  <p  ,  the 
P  Pq 

maximum  flow  rate  (maximum  <pp)  occurs  near  the  forward  wall  of  the 

vaneless  diffuser,  while  for  large  q>  it  occurs  at  the  rear  wall. 

p0 

Thus,  with  increasing  gas  flow  through  the  wheel  there  is  a  displace¬ 
ment  of  the  flow  maximum  from  the  forward  to  the  rear  wall.  In  addi¬ 
tion,  the  nature  of  the  curves  cpp  =  f(x)  attests  to  the  fact  that 
near  the  forward  wall  of  the  diffuser  the  gas  flow  across  the  section, 
for  the  majority  of  regimes,  changes  more  sharply  than  at  the  rear 
wall.  For  low  flow  rates,  reciprocal  flows  are  noted  at  the  walls. 

The  redistribution  of  the  flow  in  the  meridional  cross  section 

of  the  channel  with  increasing  ©  can  evidently  be  explained  by  the 

p0 

gradual  development  of  separation  behind  the  bend  at  the  inlet  to  the 
inner  surfaoe  of  the  oovering  disk  due  to  increased  diffusivity  in 
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this  zone  with  increasing  c  .* 

0 

8.  The  nature  of  the  curves  ol  »  f(x)  shown  in  Pig.  4  is  basically 
the  same  as  that  of  the  curves  <pp  =  f(x). 

9.  The  curve  q>  »  f  (X),  for  average  and  high  flow  rates,  qualita¬ 
tively  has  a  shape  which  is  approximately  the  opposite  of  that  of  the 
curve  <p  =  f(x),  which  also  follows  from  Formula  (10).  For  low  flow 
rates,  however,  this  connection  is  not  noted,  which  attests  tp  the 

inapplicability  of  the  Stodola  formula  when  cp  <  <p  . 

P0  P0cr 

10.  The  curves  ifP  =  f(x)  have  basically  the  same  nature  as  the 
curves  q>  -  f(X). 


11.  The  drop  in  temperature  AT*  characterizes,  for  all  practical 
purposes,  the  head  hi  and,  consequently,  to  a  certain  extent  (without 
taking  into  account  disk  and  slit  energy  losses)  the  coefficient  <p. 

The  nature  of  the  change  in  the  curves  AT*  =  f(X)  and  <p  =  f(X)  , 
should  therefore  be  approximately  analogous,  which  is  verified  by 
the  derived  experimental  curves.  Maximum  deviation  between  the  nature 
of  these  curves  is  noted  for  low  flow  rates  due  to  the  relative  increase 
in  the  leakage  of  gas  at  the  covering  disk  and  disk  losses,  which  also 
follows  from  the  theoretical  relation 

Ar*  •  -  <P 4  <’  +  P"- "  fW  =  «P(»  +  +  (W.  (19) 

.  LP 


p 

where  c  =  (A/c  )(u2/g)  =  const. 

r 

12.  The  nature  of  the  change  in  the  numbers  ,  Mjj  and  M  , 

11  ' 

and  also  M  ,  a3  a  function  of  the  wheel  regime  is  shown  in  Fig.  5. 
c0 

These  M  numbers  were  determined  as  follows  (Fig.  6): 


*  It  is  evident  that  the  flow  structure  behind  the  wheel  will 
also  be  a  function  of  the  shape  of  the  meridional  cross  section  of 
the  wheel.  We  can  also  assume  that  the  nature  of  the  change  in  the 
curves  <p  =  f(x),  a  =  f(x),  cp  =  f(x),  =  f(x)  and  AT*  =  f(x)  can  be 

influenced  by  the  leakage  near  the  covering  disk  and  through  the 
diaphragm  packing;  however,  this  question  requires  special  study. 
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Mf,  i 


w, 


ywr,  ’ 


VkgRT ,  VkgRT,  ■ 

M'  *■  -  ro'  -  M  Zs-.—L— 

-  “  / w;  t*  '■’*  T,sin|J8* 


(20) 

(21) 

(22) 


where 


Ul, 


W, 


sin  p, 
fi  *>n  t%, 


-  Vl ?  -LW?  iinPl  = 

"  K  '  ’  .  t,  sin  ft, 


M  = 


V  kgRT  ' 


fc.  ‘  T**infe, 


(23) 

(24) 


With  increased  gas  flow,  1.  e.,  <p  ,  the  numbers  M 


V 


and 


'0  ''O  1 

particularly  Increase.  This  also  follows  from  an  analysis  of 

Formulas  (20)-(23). 


The  number  Mc  depends  little  on  the  flow  rate. 


In  an  optimal 


regime  M„  -  0.17,  M  =  0.37,  M  =  0.435,  -  °*51«  In  the  regime 

cO  W1  c  W1 

i  -  0,  M  «  0.22,  M  .  «  0.4,  M  =  0.415  and  M”  =  0.66. 

J.  Cq  c 

At  the  inlet  to  the  Impeller  cascade  with  low  flow  rates  (i^  >  0) 
there  is  diffusivity  (Fig.  6),  i.e.,  /I ^  <  1.  With  increasing 

flow  rate  the  diffusivity  decreases  to  zero  and  then,  especially  at 
i^  <  0,  effusiyity  rapidly  increases,  i.e.,  /t^  >  1. 

The  influence  of  the  flow  rate  and  the  geometric  parameters  of 
the  wheel  on  the  velocity  w"  and  the  number  M"  is  evident  from 

wi 

Formulas  (23)  and  (22).  The  ratio  sin  6,/t.  sin  B  defined  by  the' angle 
of  attack  i^,  since  i^  -  f(P1,  ),  and  by  the  restriction  of  the 


*  The  superscript  ( * )  for  the  parameter  indicates  that  the 
given  parameter  is  selected  taking  into  account  restriction  of  the 
flow  by  the  vanes,  while  the  superscript  (,f)  indicates,  in  addition, 
that  aooount  is  taken  of  the  bending  of  the  flow  at  the  inlet  to  an 
angle  .  The  value  w£  is  calculated  somewhat  more  simply,  without 

taking  into  account  the  change  in  the  geometric  parameters  from  cross 
section  B-B  to  section  A- A  (see  Fig.  6).  In  the  formulas  it  has  been 

TH  and  c_  =  c. 


assumed  that  Tq  *»  T^ 
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flow  at  the  inlet  t^,  characterizes  the  diffusivity 

(sin  sin  <  l)  or  effusivity  (sin  sin  >  i)  on  entry 

into  the  inter- vane  channels  of  the  wheel.  With  increasing  gas  flow 

rate  and,  consequently,  increasing  velocities  c  and  c  ,  the  numbers 

z0  rl 

,  Mr  will  inorease,  as  can  be  seen  from  Formulas  (20),  (21), 

and  (22).  But  the  intensity  of  the  increase  in  ,  compared  with 

M  ,  in  this  case  will  be  greater,  since  the  ratio  sin  B./sin  0  will 
wi  ‘  1  gl 

also  increase  due  to  an  increase  in  the  angle  0^  Therefore  with 

increasing  flow  rate  the  diffusivity  at  the  inlet  to  the  inter-vane 

channels  should  decrease.  Theoretically,  it  follows  from  Formulas 

(21),  (22),  and  (23)  that  with  abnormally  high  restriction  (very  low 

T^),  when  the  ratio  sin  0^/^  sin  Pg  >  1  and,  consequently,  , 

the  diffuser  zone  will  generally  be  absent  at  the  inlet,  even  with 

large  positive  angles  of  attack. 

The  point  of  transition  from  a  diffuser  zone  to  an  effuBer  one 

is  determined  by  the  condition 


t,  sin  B 


*1 


or,  after  transformation. 


9p.M  — 


t,  sin  p 


V\- 


^iVi  A 

^«Yo  '  A  ' 


(25) 


.(26) 


where  9  M  is  the  value  of  the  flow  factor  9  for  which  w. 


1* 


'P0M  'Po 

The  point  of  intersection  of  the  curves  1^  and  (Fig.  5) 

strictly  speaking,  is  determined  by  the  condition  or 


f  T, 


s*n  Pi _  .  /  'J _ 1 

T,  »in  p  \  r\ 


(27) 


The  ratio  when  is  approximately  equal  to  unity, 

which  was  assumed  to  be  the  case  during  the  calculations.  Therefore 
the  point  of  intersection  of  the  curves  and  ^  is  treated  as  the 
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point  of  transition  from  a  diffuser  zone  to  an  ef fuser  one.  This 
transition  point  is  always  in  the  region  i^  >  0,  since  <  1*  and, 

consequently,  w^  =  w^  only  when  sin  p^/sin  Pg  <  1,  i.e.,  <  Pg 

and  1^  >  0.  The  transition  point  can  be  found  when  i1  =  0  only 
theoretically  —  in  the  absence  of  restriction,  when  «  1. 

Conclusions 


1.  Investigation  has  made  it  possible  to  obtain  a  more  complete 
understanding  of  the  efficiency  of  operation  of  the  wheel  of  a  centri¬ 
fugal  compressor  and  to  quantitatively  associate  individual  experi¬ 
mental  data  with  the  '"esults  of  theoretical  analysis. 

2.  The  attained  efficiency  =  0.9  is  maximum  for  wheels  with 

geometric  parameters  similar  to  the  investigated  wheels.  For  other 

geometric  parameters,  in  particular  the  angles  p  and  B  ,  other 

g2  gl 

peripheral  velocities  Ug,  etc.,  the  efficiency  level  may  change.  For 
example,  the  transition  to  a  two-stage  cascade,  as  experimental  data 
have  shown  (see  p.  44),  makes  it  possible  to  raise  the  efficiency  of 
the  wheel.  A  similar  statement  relative  to  the  dependence  of  the 
obtained  experimental  data  (mainly  quantitative)  on  the  initial  geo¬ 
metry  of  the  investigated  wheels  may  be  made  for  the  other  criterional 
parameters. 

J.  The  experiments  were  conducted  using  wheels  made  according 
to  the  generally-accepted  recommendations  for  -typical  wheels  of  sta¬ 
tionary  centrifugal  compressors  of  limited  performance.  Therefore 
the  obtained  data  on  the  values  of  the  various  criterional  parameters 
of  the  wheel  can  be  used  direotly  mainly  in  designing  the  indicated 
compressors. 
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Fig.  2.  Diagram  of  the  blading 
of  the  ETaK-1  stand. 
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Fig.  Characteristics  of 
the  wheel  z^  =  Zg  =  16: 


^  T'  X" 


pm 
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l8. 


't 

i 


Fig.  5.  Characteristics 
of  the  wheel  according 
to  number  MCCz^  =  zg  =  1 6). 


Fig.  6.  Inlet  velocity 
triangles  of  a  centrifugal 
wheel  with  and  without 
taking  into  account  the 
restriction  and  bending" 
of  the  flow. 
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THE  INFLUENCE  OP  THE  NUMBER  OF  VANES  ON  THE  EFFICIENCY 
OF  THE  CENTRIFUGAL  WHEEL  WITH  A  SINGLE-STAGE  CASCADE 


S.  A.  Anisimov,  F.  S.  Rekstin,  and  K.  P.  Seleznev 

The  operating  efficiency  of  a  wheel  depends  on  the  type,  the 
relative  dimensions,  and  the  operating  regime  of  the  wheel.  In  addi¬ 
tion  to  the  M  and  Re  numbers  and  the  angle  of  attack  i.  (or  cp  )  for 

x  Po 

a  given  medium  with  a  specific  heat  ratio  k,  the  basic  parameters 
which  define  the  operation  of  the  wheel  are  the  geometric  character¬ 
istics  pg  ,  b2/D2,  z2,  Pg^,  D1/D2,  Dq/Dq,  F/Fq,  the  shape  of  the 

vane  profile,  and  others.* 

In  this  article  we  analyze  the  results  of  an  experimental  study 
of  the  influence  of  the  number  of  vanes  of  a  centrifugal  wheel  on  its 
characteristics  and  on  the  flow  structure  behind  it.  The  work  was 
conducted  on  the  ETsK-1  test  stand  at  the  Compressor  Machine  Labora¬ 
tory  of  the  Kalinin  Leningrad  Polytechnic  Institute. 

The  diagram  of  the  blading  of  the  stand  is  shown  in  Fig.  2  on 

p.  16. 

The  tests  were  conducted  using  ten  shrouded  wheels  with  entirely 


*  The  arbitrary  designations  used  in  this  article  are,  with  very 
few  exceptions,  the  generally-accepted  ones. 
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milled  vanes  having  an  angle  0  =  49° •  The  wheels  differed  only  in 

82 

the  number  of  vanes:  z  =  z^  =  z2  =  4,  6,  8,  10,  12,  14,  16,  20,  24, 
and  28.  The  basic  geometric  parameters  of  these  wheels  have  been 
given  in  a  previous  article  [l].  The  diffusivity  of  the  channels  and 
the  wheel  F^/Fj,  FJ^/F^  and  F^/F^,  the  relative  pitch  of  the  equivalent 
straight  cascade  t/l,  the  coefficients  of  flow  restriction  by  the 
vanes,  and  the  angles  of  exposure  of  the  inter-vane  channels  in  the 
plane  of  rotation  at  the  inlet  x^  and  v^,  and  at  the  outlet  x2  and  v2, 
for  a  series  of  test  wheels  are  given  in  Fig.  1.  Here, 

•  ji  0,6,  — *,4,- 


T*  =  “  = 

K 


a, 

sin  B 

_ -8.!.-  —  1 


T.  = 


2 

77' 


nD1bl 
.  nDtb,  —  z.6,  - 


Mi 


nDtb. 


=  i 


sin 

nO,  sin 


z  In 


0. 

0. 


^2  _  ^  «TS  *‘n  ftg,  _  ^2  >ln  , 
F\  >■' if,  sinfc,  ~  '  sin  ft,/ 


=  ( 

■■  \ 


COS 


SCO/,  0.  -  \ 

=  V1  — ^cosfe.) 


(1) 

(2) 

(5) 

(4) 

(5) 

(6) 


where 


/?,  =r 


1  2(«.  cos  fe,  — KiCOsfc,.)' 


From  the  graphs  in  Fig.  1  it  follows  that  when  z  =  4-28  the 
coefficients  x^  =  0.91-0. 39  and  t2  =  0.97-0.78,  (when  z  =  16,  x^  =  0.66 
and  x2  =  O.87).  As  can  be  seen  from  these  numbers,  the  investigated 
wheels  had  great  restriction,  which  is  characteristic  of  limited-per¬ 
formance  centrifugal  wheels.  For  comparison,  we  give  the  values  of 
x  for  the  same  type  of  wheel  on  two  machines  of  the  Neva  Machine- 
Construction  Factory:  the  K-100-61-1;  x^  **  O.73,  x2  «  O.87,  and 
K-500-61-1;  t±  *»  O.73-O.7 6,  x2  -  O.89. 
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1.  The  Problem 


In  solving  the  problem  of  selection  of  the  best  number  of  vanes 
for  the  wheel  we  are  usually  guided  by  the  following  concepts;.  As  is 
known  [5,  4,  5],  the  theoretical  head  h^  whioh  the  gas  attains  In  the 


wheel  channels  when  c..  =  0  is  determined  from  the  formula 


u. 


hT  — 


(7) 


where  <p2  =  cu2/u2* 

Using  as  our  basis  the  Stodola  formula,  which  is  often  used  in 
the  theory  of  stationary  centrifugal  compressors  to  calculate  the 
influence  of  finite  number  of  vanes,  we  get 


<pa  =  1  _  JL  Sin  ~ 


(8) 


From  this  it  follows  that  with  increasing  z2  there  is  an  increase 
in  q>2  and,  consequently,  in  h,p.  In  deriving  Formula  (8)  we  took  into 
account  only  the  effect  of  the  relative  eddy  on  the  velocity  value  c 
because  of  the  finite  number  of  vanes. 

However,  with  the  flow  of  a  real  gas  in  a  centrifugal  wheel  with 
a  finite  number  of  vanes,  in  addition  to  the  action  of  the  relative 
eddy,  there  are  a  number  of  other  factors  which  influence  the  velocity 
triangle  and,  consequently,  the  value  of  <p2* 

With  a  considerable  reduction  in  z  the  load  on  each  vane  increases, 
the  pressure  and  velocity  gradients  are  raised  on  the  surface  of  the 
vane,  and  the  intensity  of  the  diffuser  zone  on  the  profile  increases. 
There  is  a  simultaneous  increase  in  the  angle  of  exposure  of  the 
inter-vane  channel  and  local  values  of  the  Mach  number.  These  all 
result  in  the  appearance  and  development  of  separation  zones  which 
disrupt  flow  in  the  channel,  cause  a  decrease  in  the  head  due  to 
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decreased  lifting  force  and,  as  a  result,  additional  deviations  in 
the  flow  at  the  outlet,  and  also  lead  to  decreased  efficiency.  The 
separation  zones  restrict  the  outlet  cross  section,  increase  the  flow 
velocity,  and,  consequently,  decrease  the  coefficient  <p2  according  to 
Formula  (8).  From  what  has  been  said  it  is  clear  that  with  a  small 
number  of  blades  the  head  will  decrease  more  sharply  than  follows 
from  theoretical  relations  (7)  and  (8). 

On  the  other  hand,  with  an  increased  number  of  blades  their  sur¬ 
face  greatly  increases,  i.e.,  frictional  losses  increase,  which  results 
in  a  drop  in  the  useful  head  and  the  efficiency  of  the  wheel.  There 
is  simultaneously  increased  diffusivity  of  the  wheel  channel  (see 
Fig.  l)  and  a  great  increase  in  the  gas  velocity  at  the  inlet  to  the 
channel  due  to  sharp  restriction  of  the  flow  (small  values).  These 
factors  also  cause  the  development  of  separation  zones  and,  as  a 
result,  a  decrease  in  the  theoretical  and  useful  heads. 

From  what  has  been  said  it  is  clear  that  the  head  and  the  effi¬ 
ciency  of  the  wheel  should  be  maximum  for  some  average  number  of  vanes. 

Another  important  parameter  which  determines  the  efficiency  of 
the  wheel  is  the  width  of  the  operating  zone  of  the  characteristic. 

Let  us  examine  this  question  for  the  case  of  the  flow  of  an 
ideal  gas.  Theoretical  investigations  have  shown  [3,  4]  that  under 
specific  conditions  (an  lnviscid  gas,  unstressed  entry,  low  output) 
the  gas  velocity  w  on  the  working  side  of  the  blade  (the  side  facing 
the  flow),  because  of  the  Influence  of  the  relative  eddy,  can  drop 
to  zero  or  even  become  negative.  This  latter  will  indicate  the  appear- 
anoe  of  invars®  flows  on  the  acting  surface  of  the  vane  and,  conse¬ 
quently,  unstable  gas  motion.  Such  flow  motion  causes  the  appearance 


acting  surface  of  the  blade  decreases  to  zero  determines  the  boundary 
of  stable  wheel  operation  and  is  called  the  critical  output.  This 
output  is  associated  with  the  optimum  output  by  the  following 

relation  [4]: 


2.nt  sin*fe. 


(9) 


where  a  is  the  width  of  the  vane  channel  in  the  wheel  revolution  plane. 

Prom  this  formula  it  follows,  in  particular,  that  with  an  Increased 
number  of  vanes  Vgr/Vgpt  decreases,  i.e.,  the  acting  zone  of  the 
characteristic  Increases.  However,  the  flow  structure  of  an  ideal 
gas  in  a  wheel  in  a  design  and  nondesign  regime,  as  experimental 
studies  have  shown,  differs  considerably  from  the  theoretical  structure 
[3,  4].  This  is  manifested,  in  particular,  by  the  fact  that  inverse 
flow  does  not  occur  along  the  acting  side  of  the  vane,  even  with  very 
low  outputs.  And,  conversely,  on  the  nonacting  side  of  the  vane  there 
is  always  breakaway,  which  is  the  more  noticeable  the  lower  the  out¬ 
put.  These  facts  make  it  necessary  to  use  with  care  those  conclusions 
from  Formula  (9)  dealing  with  the  influence  of  the  number  of  vanes 
on  the  value  of  the  acting  zone  of  the  characteristic.  The  difference 
between  the  actual  flow  pattern  and  that  of  an  ideal  gas  is  explained 
by  the  appearance  of  separation  in  the  boundary  layer  on  the  diffuser 
segments  of  the  nonacting  surface  of  the  vanes  and  by  their  develop¬ 
ment,  which  results  in  increased  average  relative  velocity  w. 

Finally,  an  increase  in  z_  obviously  complicates  the  design  and 
construction  of  the  wheels,  i.e.,  it  increases  their  cost. 

The  relation  between  the  actual  flow  structure  in  the  wheel 
channels  and  the  number  of  vanes,  as  is  clear  from  the  preceding,  is 
extremely  complex  and,  at  least  at  present,  has  not  been  precisely 
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determined  theoretically.  Thus  the  problem  of  selecting  the  number 
of  vanes  should  be  solved  by  taking  into  account  theoretical  and 
experimental  data.  For  practical  purposes,  the  number  of  vanes  should 
be  selected  depending  on  which  wheel  characteristics  are  used  as  the 
determining  ones,  viz.,  efficiency,  head,  width  of  acting  zone  of  the 
characteristic,  flow  structure  behind  the  wheel,  reliability,  feasi¬ 
bility,  etc.  This  latter  results  from  the  fact  that  at  present  it 
is  actually  impossible  to  create  a  wheel  to  satisfy  simultaneously 
all  these  various  requirements.  In  other  words,  a  wheel  may  be 
optimum  only  with  respect  to  certain  characteristics.  This  is  the 
difficulty  in  selecting  the  number  of  vanes  and/  in  general,  the 
arbitrary  nature  of  this  concept. 

In  stationary  centrifugal  compressors  in  the  past,  the  wheels 
were  made  with  a  relatively  large  number  of  vanes  (z  ^  28).  At 
present,  fewer  vanes  are  used.  For  example,  at  the  Neva  Machine-Con¬ 
struction  Factory,  for  typical  "compressor"  (B  =  45°  and  c  /u0  = 

g2  r2 

-=  0.24)  and  "pump"  (B  =  20°  and  c  /u,-  =  0.12)  wheels,  18  and  8  vanes, 
respectively,  are  used  [6]. 

Various  analytical  dependences  are  used  In  the  literature  to 
select  the  appropriate  number  of  vanes.  These  dependences,  and  also 
the  results  of  using  them  to  calculate  the  recommended  number  of 
vanes  for  the  wheels  we  investigated,  are  given  in  the  table.  As  can 
be  seen  from  the  table,  the  recommended  number  of  vanes  varies  within 
a  rather  broad  range  —  from  8  to  22  —  depending  on  which  formulas 
were  used  in  the  calculations. 

The  experimental  data  in  Soviet  and  foreign  literature  on  the 
Influence  of  the  number  of  vanes  of  a  centrifugal  wheel  on  its  charac¬ 
teristics  are  incomplete  and  pertain  to  studies  of  wheels  which  are 
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not  typical  for  stationary  centrifugal  compressors  [7,  14],  or  to 
determination  of  the  total  characteristics  of  centrifugal  stages  with 
wheels  having  various  z  [153.  The  literature  contains  no  detailed 
data  on  the  influence  of  the  number  of  vanes  on  the  flow  structure 
behind  the  impeller. 

2.  Research  Methods 


The  method  of  investigating  the  influence  of  the  number  of  vanes 
of  a  wheel  on  the  efficiency  of  its  operation  is  based  mainly  on  the 
method  of  investigating  the  operation  of  individual  wheels  given  in 
previous  papers  [l,  2],  The  main  difference  is  that  if  when  investi¬ 
gating  individual  wheels,  as  an  independent  parameter  we  use  <p  or  X, 

p0 

then  when  investigating  the  influence  of  the  number  of  vanes  it  is 
natural  to  use,  as  the  independent  parameter,  the  number  of  vanes  and 
also  the  coefficient  t^,  the  relative  pitch  t/l,  and  the  angle  of 
exposure  of  the  vane  channel  v^,  uniquely  associated  with  z  for  the 
investigated  wheels  and  giving  a  more  complete  idea  of  the  conditions 
of  flow  entry  into  the  wheel  channels.  In  addition  to  the  previously 
used  criterlonal  parameters  (see  the  footnote  on  p.  J>)  r)*^,  ,  <P,  ft, 

p,  kfi,  kQ,  etc.,  to  investigate  the  influence  of  the  number  of  vanes 
on  the  width  of  the  acting  zone  of  the  characteristic  of  the  wheel 
we  introduced  the  auxiliary  criterlonal  parameter,  the  coefficient  k^.: 


yp.opt  Vp.  cr 
'h’.or 


(10) 


where  <p  and  cp, 
^0 


p^cr  pQopt 


are  the  gas  flow  factors  <p  for  tJr*  and 

p  q  max 


max'  respectively. 


When  ky  increases  there  will  be  expansion  of  the  acting  zone  of 
the  wheel  characteristic. 
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3.  Investigation  Results 


1.  The  curve  max,  shown  In  Fig.  2,  is  smooth  for  a  broad 
range  of  z  (z  ■  8-20). 

When  z  <  8  there  Is  a  noticeable  droj}  In  the  efficiency  due  to 
deterioration  of  the  flow  structure  in  the  wheel  channels,  associated 
with  separation  phenomena.  With  such  a  small  number  of  vanes  the  load 
evidently  increases  greatly  on  the  vane,  and  the  exposure  angle  of 
the  vane  channel  increases  considerably,  which  results  in  the  vigorous 
development  of  separation  phenomena.  In  this  case,  in  the  vane  chan¬ 
nels  there  occur  local  pressure  and  velocity  gradients. 

When  z  >  20  there  is  a  certain  decrease  in  the  efficiency  due  to 
the  strong  restriction  of  the  flow  at  the  inlet  (the  velocity  w^  and 

the  number  M1'  ) ,  an  increase  in  the  diffusivity  of  the  channels  and 
1  , 

the  wheel  (F£/F^,  F2/F£,  and  Fg/F^  increase),  and  the  influence  of 
friction. 

The  maximum  is  in  the  region  z  =  10-18.  In  this  case  the 

angles  of  attack  are  within  the  limits  i1  «  3-7°. 

The  numerical  value  of  the  angle  of  attack  depends  on  the  shape 

of  the  leading  edge  of  the  profile.  The  angle  of  attack  is  calculated 

from  the  angle  which  is  determined  by  the  direction  of  the  camber 
gl 

line  at  radius  R1,  while  the  leading  edge  of  the  blade  is  practically 
symmetrical.  Depending  on  the  shape  of  the  leading  edge,  optimum 
flow  can  occur  for  several  values  of  i^.  For  example,  with  tapering 
of  the  leading  edge  of  comparatively  thick  completely-milled  vanes 
from  the  concave  side,  unstressed  flow  will  be  attained  at  i^  >  0. 

In  addition,  we  should  bear  in  mind  the  influence  on  the  true 
value  of  the  angle  of  attack  of  leakages  near  the  covering  disk, 
although  their  influence  is  quite  insignificant  and  is,  on  the  average. 
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for  the  examined  wheels  (radial  gap  of  0.25-0.35  mm)  in  optimal 
regimes  Ai1  £  30-40*.  The  leakages  decrease  the  algebraic  value  of 
the  angle  of  attack. 

The  above  phenomena  have  remained  practically  unstudied;  there¬ 
fore  it  is  very  difficult  to  estimate  the  actual  extent  of  the  influ¬ 
ence  of  these  factors.  Special  experimental  and  theoretical  operations 
are  required  to  solve  these  problems. 

2.  The  shape  of  the  curves  „  =  f(z)  and  fflD„  =  f(z)  shown 

in  Fig.  2  is  similar  to  that  of  the  curve  r[*^  max  =  f(z).  The  decrease 

in  the  value  of  and  q>  „  with  small  and  large  z  can  be  explained 

max  max 

by  the  same  reasons  as  the  decrease  in  Ty*^  max.  The  maximum  and  <p 
are  in  the  region  z  =  14-20. 

3.  Figure  3  gives  individual  characteristics  of  wheels  with 
various  numbers  of  vanes.  From  the  graphs  it  follows  that  since  the 
curves  and  ip*  are  more  sloping  with  decreasing  number  of  vanes, 
consequently  with  decreasing  z  the  zones  of  economical  and  stable 
wheel  operation  increase. 

The  increased  zones  of  economical  and  stable  operation  with 
decreasing  number  of  vanes  can  evidently  be  explained  as  follows. 

When  i^  *»  0  flow  past  the  leading  edges  is  not  characteristic;  the 
wheel  characteristics  are  determined  by  conditions  of  gas  flow  in  the 
wheel  channel.  When  i^  ^  0  and  their  absolute  values  are  reasonably 
high,  the  maximum  local  M  numbers  and  maximum  local  diffusivity  occur 
with  gas  flow  near  the  leading  edges.  The  nature  of  the  flow  past 
these  edges  for  any  z  should  evidently  be  similar,  since  for  all  the 
investigated  wheels  their  shape  is  identical,  while  the  value  of 
maximum  diffusivity  and  the  maximum  M  number  in  this  region  should 
be  determined  mainly  by  the  average  value  of  the  M  number  near  the 
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leading  edges.  From  the  graph  in  Fig.  6  we  see  that  with  an  increase 

in  z  the  number  M"  continually  increases,  i.e.,  with  an  increase  in 

z  breakaways  begin  to  develop  at  lower  absolute  values  of  the  singles 

of  attack  than  we  observed  in  the  investigated  wheels. 

With  an  increase  in  z  the  efficiency  maxima  shift  toward  smaller 

<p  .  while  the  head  coefficient  maxima,  on  the  other  hand,  shift 
p0 

toward  higher  cp  ,  i.e.,  in  this  case  we  observe  a  certain  decrease 
p0 

in  cp  „  .  and  an  increase  in  CP  This  indicates  that  with  an 

Pq°Pu  Po  1 

increasing  number  of  vanes  the  width  of  the  acting  zone  of  the  wheel 
characteristic  defined  by  the  coefficient  ky  (Fig.  2)  decreases.  This 
latter  does  not  confirm  theoretical  formula  (9). 

The  shift  in  the  maxima  of  TtfU  and  V'*  vs.  z  can  evidently  be 
explained  by  the  influence  exerted  by  the  change  in  the  value  of  flow 
constriction  at  the  inlet  on  the  ratio  between  and  cr. 

The  critical  number  cr  is  the  value  of  Mw  at  which  at  any 
point  on  the  surface  of  the  profile  the  gas  velocity  attains  the  local 
velocity  of  sound.  The  value  of  cr  depends  basically  on  the  geo¬ 
metry  of  the  cascade  design  and  the  angle  of  attack.  To  assure 
normal  operation  of  the  cascade  we  must  have  the  inequality  ^  'V^cr' 
Otherwise,  on  the  vane  surface  there  arises  a  supersonic  zone  and, 
as  a  result,  there  are  compression  shocks  which  cause  the  vigorous 
development  of  breakaway  phenomena  and  a  sharp  deterioration  of  the 
flow  aerodynamics  in  the  wheel  channels. 

The  value  of  the  local  velocity  in  the  vicinity  of  the  leading 
edge  depends  mainly  on  the  shape  of  the  edge,  the  angle  of  attack, 
and  the  value  of  the  velocity  ahead  of  the  cascade  and  in  the  throat 
of  the  cascade  channel.  The  local  velocity  for  the  given  geometry 
and  specific  velocity  ahead  of  the  cascade  is  minimum  when  i^  «  0 
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and  increases  when  i„  ^  0.  Por  the  given  angle  of  attack  the  local 
velocity  increases  with  increasing  constriction,  i.e.,  with  increasing 
z  (except  for  the  case  of  very  small  z_)  due  to  increased  velocity  in 
the  narrow  channel  cross  section. 

There  are  practically  no  data  on  the  value  of  cr  for  circular 
cascades.  In  this  article  let  us  give  .an  analysis,  using  the  concept 
of  the  critical  number  cr,  by  which  we  will  understand  that 
number  [l]  for  which  somewhere  on  the  profile  there  occurs  a  local 
velocity  equal  to  the  local  velocity  of  sound.*  The  number  cr  has 
a  maximum  value  at  i^  «  0  and  decreases  when  i1  >  0  and  i1  <  0. 

Regimes  corresponding  to  max  lie  in  the  region  of  small 
angles  of  attack  close  to  zero.  In  this  case  there  are  optimum  con¬ 
ditions  for  vane  profile  flow,  due  to  which  the  number  cr  has  a 
high  value  and  changes  but  little  with  slight  changes  in  the  angle 

of  attack,  i.e.,  in  this  case  preservation  of  the  Inequality  M"  < 

W1 

<  cr  will  depend  mainly  on  the  change  in  the  value  of  and, 
consequently,  with  constant  cascade  geometry,  only  on  the  value  of 
the  coefficient  <p. 


p0 

Let  us  write  the  expression 


M.,= 


T 


sin  Pi 
Ti*in  f 


(/  unm] 


\fkgm\ 


t,sin  \f kgRT] 

From  thi3  it  follows  that  with  Increasing  velocity  c 


(11) 

and. 


consequently,  <p  ,  and  with  Increasing  flow  restriction  (i.e.,  with 
p0  , 
decreasing  ),  the  velocity  w"  and  the  number  M"  will  increase. 

X  X 


*  The  critical  number  1^  cr  should  be  considered  arbitrary, 

since  for  not  very  high  Ug  and  small  or  average  values  of  _z  It  is 

possible  that  local  supersonic  velocity  is  not  attained.  In  this 
case,  cr  characterizes  the  regime  for  which  due  to  high  local 

diffusivity  breakaways  vigorously  develop. 
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Thus  the  influence  of  the  flow  rate  and  flow  restriction  on  M"  is 

W1 

qualitatively  identical.  Or,  in  other  words,  in  order  that  the 

number  Bl"  not  change  with  increasing  restriction  due  to  an  increase 
W1 

in  the  number  of  vanes  it  is  necessary  that  the  gas  flow  rate  (i.e., 

<P_  )  decrease.  Therefore,  in  the  regimes  r )*n  ma.  when  preservation 
Pq  3,0  max 

of  the  inequality  <  M”  cr  depends  mainly  on  the  value  of  Mm  , 

with  increasing  _z  there  must  be  decreased  gas  flow  to  prevent  an 

increase  in  .  From  this  it  is  understandable  that  with  increasing 

z  the  coefficient  ©  .  should  decrease. 

Pq°P 

Regimes  corresponding  to  lie  in  the  region  of  large  positive 

angles  of  attack  and  have  correspondingly  poor  conditions  of  vane 

profile  flow,  due  to  which  the  number  cr  has  a  relatively  low  value 

and  depends  greatly  on  the  angle  of  attack,  i.e.,  in  this  case  with 

constant  cascade  geometry  preservation  of  the  inequality  M"  <  ^  cr 

will  be  determined  mainly  by  the  value  cr,  which  changes  abruptly 

even  with  insigifnicant  changes  in  the  angles  of  attack  or  cp  . 

p0 

Increased  flow  rate  in  this  case  leads  to  a  decrease  in  i^  and  there¬ 
fore  causes  a  considerable  increase  in  CI>.  In  this  case  there  Is 
also  a  certain  increase  in  ;  however,  we  must  assume  that  an  increase 

in  the  flow  factor  influences  an  Increase  In  the  number  M"  to  a  much 

W1 

lesser  extent  than  an  increase  In  cr,  i.e.,  the  inequality  < 

<  B^  cp  is  retained.  With  decreasing  flow  rate,  i.e.,  an  increase 
in  the  angles  of  attack,  the  opposite  should  be  observed.  In  this 
event  the  Inequality  M”  <  B4"  is  not  preserved. 


Thus,  in  regimes  when  the  decisive  factor  is  BT  cr,  reten¬ 


tion  of 


the  inequality  <  PT  cr  with  increasing  z  (and,  consequently. 


the  number  M’’  )  is  possible  only  when  the  flow  rate  changes  such  that 
W1 


;he  number 


\or 


increases.  This  will  occur  with  a  decrease  in  the 
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angle  of  attack,  l.e.,  with  an  increase  in  cp  [1].  From  this  we 

p0 

see  that  with  increasing  z  the  coefficient  <p  corresponding  to  . 

—  Pq  max 

l.e.,  <p  cr,  should  Increase,  which  is  confirmed  by  experimental  data 
obtained  (see  Figs.  2  and  3). 

In  addition,  from  Formula  (ll)  it  follows  that  the  qualitative 

influence  of  on  the  number  M"  when  i.  =  const,  i.e.,  when 
«*•  w1  ± 

sin  p./sinp  =  const,  is  opposite  to  the  influence  of  t. ;  therefore 

we  can  assume  that  for  wheels  similar  to  those  investigated  an  increase 

in  the  velocity  u^,  when  its  value  is  reasonably  high,  should  cause 

qualitatively  the  same  type  of  shift  of  the  maxima  of  q*^  and  f*  as 

a  decrease  in  t^,  i.e.,  an  Increase  in  the  number  of  vanes  z.  In 

other  words,  to  expand  the  acting  zone  of  the  wheel  characteristic 

we  must  decrease  the  peripheral  velocity  of  the  wheel  u^.  This  latter 

conclusion  naturally  requires  experimental  verification. 

A  change  in  the  angle  £  ,  when  the  remaining  geometric  para- 

gl 

meters  of  the  wheel  are  retained,  results  in  a  change  in  the  optimum 

flow  rate,  l.e.,  the  optimum  value  of  the  coefficient  <p  .  Using 

p0 

wheels  with  large  p  it  is  evidently  possible  to  displace  the  zone 

S1 

of  maximum  efficiency  to  the  region  of  increased  flow  rates. 

The  results  obtained  on  the  shift  in  the  maximum  of  the  curves 

f*  =  f(cp  )  with  a  change  in  z  agree  with  the  experimental  data  of 
p0 

Kearton  [14]  and  Ris  [15].  The  shift  in  the  maximum  of  the  curves 

q**  =  f(q>  )does  not  agree  with  the  experimental  data  of  these  authors; 
ao  p0 

thi3  latter  can  evidently  be  explained  by  the  influence  of  the  fixed 

elements,  since  our  data  pertain  to  a  wheel,  while  the  other  data 

[14,  15]  correspond  to  a  stage  and  a  unit,  respectively. 

4.  The  coefficients  of  relative  flow  nonuniformity  with  respect 

to  the  internal  k^  and  adiabatic  k£  heads  (Figs.  4  and  5)  are 

a  ad 
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practically  independent  of  the  number  of  vanes,  both  for  t)*^  and 
il  "  °>  “d  for  Kax‘ 

The  relative  flow  nonuniformity  with  respect  to  flow  rate  and 
angle,  characterized  respectively  by  the  coefficients  kQ  and  ka  (Figs. 
4  and  5)>  depends  on  z,  particularly  for  ^ax*  This  is  associated 
with  the  fact  that  the  operational  regimes  corresponding  to  lie 

close  to  the  surge  zone  where  breakaway  phenomena  are  particularly 
greatly  manifested. 

For  regimes  max  and  i^  =  0  a  certain  increase  in  kQ  and  kQ 
is  observed  with  a  decrease  of  z  <  8,  due  to  the  development  of  break¬ 
away  phenomena  because  of  increased  load  on  the  blades  and  increased 
velocity  and  pressure  gradients,  and  with  an  increase  in  z  >  14,  due 
to  the  increased  diffusivity  of  the  channels  and  the  wheel  and  also 
increased  friction  in  the  channels. 

The  absolute  nonuniformity  with  respect  to  the  angle  Act  =  amax  - 
-  otmln  (Figs.  4  and  5)  when  z  =  10-20  is  minimum  and  remains  practi¬ 
cally  constant,  but  increases  considerably  when  z  =  8.  In  our  test 
wheels  the  absolute  nonuniformity  with  respect  to  the  angle  was  6-10°, 
i.e.,  the  variation  of  the  angle  a  about  the  average  value  was  +(?-5)° 
even  for  the  regimes  max  and  i1  =  0.  The  presence  of  a  consider¬ 
able  absolute  nonuniformity  with  respect  to  the  angle  even  for  design 
regimes  should  be  taken  into  consideration  when  designing  vaned  dif¬ 
fusers  and  particularly  for  the  shape  of  the  leading  edges  of  the 
vanes. 

5.  The  nature  of  the  change  in  the  degree  of  reaction  ft  with  a 
change  in  z  is  determined  basically  by  the  nature  of  the  change  In 
<p  [l],  while  the  degree  of  reaction  p  is  determined,  in  addition,  by 
energy  losses  (Fig.  4). 
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The  values  of  Cl  and  p  decrease  greatly  with  an  increase  in  z  to 
6-8.  With  a  further  increase  in  z  the  decrease  in  fl  and  p  occurs 
considerably  more  slowly  (when  T)*^  max  and  i^  =  0),  while  for 
the  values  of  Cl  and  p  in  this  region  remain  practically  constant.  For 
all  numbers  of  vanes,  as  follows  from  the  theory  Cl  >  p  [l], 

6.  A  comparison  of  the  curves  9  =  f(z)  and  ®  =  f(z)  [l]*  for 


Kb 


max 


(Fig.  4)  and  also  for  1.  =  0  shows  that  they  are  qualitatively 


identical,  but  the  numerical  value  of  <p„  is  less  than  the  experimental 
value  of  <p.  With  an  increase  in  z  the  values  of  <p  and  cp  converge. 

“  C 

These  data  agree  with  those  derived  by  Kearton  [14]. 

7-  The  dependence  of  the  numbers  M  ,  M  ,  M"  and  M  on  the 

Cq  w±  w±  c 

number  of  vanes  for  T]*-v  ,  i.  =  0,  and  ip*  is  shown  in  Fig.  6.  The 

cio  mux  -L  max 

nature  of  the  curves  of  M  and  M  is  determined  by  the  change  in 

cQ  w1 

the  coefficient  ®  corresponding  to  the  regimes  Ty*-x  _a.  i.  =  0  and 
]Pq  ao  max  .l 

_  for  various  z.  The  nature  of  the  curve  for  all  regimes  is 
max  —  c 

determined  by  the  influence  of  the  number  of  vanes  on  9  and  Cl  [l]. 

The  value  of  with  increasing  z  basically  increases,  most  intensely 
in  the  region  z  =  4-10,  and  very  slowly  for  z  >  14-16. 

When  i^  «  0  the  number  continually  increases  with  increasing 
z,  since  restriction  increases  (t.  decreases).  For  ip*  the  value  of 
with  increasing  z  increases  considerably  due  to  increased  restric¬ 
tion  as  well  as  due  to  the  shift  of  to  the  region  of  high  flow 

rates,  where  effusivity  Increases  at  the  cascade  inlet. 

The  value  of  M"  for  n*^  ln  the  re6lon  z  £  16  remains  approxi- 
w^  ‘ao  max 

mately  constant.  This  constancy  of  the  value  of  with  a  change 

W1 


*  The  coefficient  9.  =  c/u0  is  found  by  converting  from  diameter 

O  U  d 

D2  to  the  diameter  of  the  control  section  D  of  the  coefficient  92  = 

-  c„  /uo  defined  by  the  Stodola  formula. 

Ug  * 
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supports  the 


in  the  number  of  vanes  in  the  region  z  £  16  for  max 

previously  expressed  concepts  on  the  influence  of  the  number  on 

the  position  of  the  maximum  efficiency.  With  a  further  increase  in 

the  number  of  vanes,  preservation  of  the  constancy  of  would  result 

in  'a  great  reduction  in  the  flow  rate  and  a  great  increase  in  the 

angle  of  attack,  in  which  case  the  number  cr  would  noticeably 

decrease.  Therefore  the  maximum  efficiency  for  z  S  16  is  attained 

for  approximately  constant  angles  of  attack  corresponding  to  the 

regime  in  which  cr  is  still  large  and  approximately  equal  to  its 

value  when  i^  =  0.  For  somewhat  larger  angles  of  attack  the  value  of 

cr  will  evidently  decrease  considerably.  It  is  characteristic 

that  when  z  >  16  m  begins  to  gradually  decrease,  which  is 

evidently  associated  with  an  increase  in  M"  ,  and  which  supports  our 

W1 

conclusions. 

The  coefficient  <p  M  (the  flow  factor  cp  ,  for  which  w.  =  w") 

Per  po  j-  i 

decreases  with  Increasing  restriction  due  to  Increasing  z,  as  can  be 
seen  from  Fig.  6.  Thus,  with  an  Increase  in  the  number  of  vanes  the 
transition  point  from  a  diffusion  zone  to  an  effusion  zone  shifts 
toward  lower  flow  rates.  Or,  In  other  words,  with  increasing  £  the 
region  of  operational  regimes  for  which  there  is  diffusion  Inlet  into 
the  vane  channel  of  the  wheel  decreases. 


Conclusion 

1.  Our  studies  have  permitted  a  qualitative  and  quantitative 
estimate  of  the  Influence  of  the  number  of  vanes  £  on  the  operational 
efficiency  of  a  wheel.  Obviously,  however,  the  specific  quantitative 
data  on  the  influence  of  the  number  of  vanes  (e.g.,  the  number  of 
vanes  corresponding  to  maximum  efficiency)  are  relative  in  nature. 
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This  is  due  to  the  fact  that  the  wheel  efficiency  depends  on  many 
factors,  and  therefore  quantitative  data,  strictly  speaking,  are  valid 
only  for  specific  design  relationships  for  a  wheel,  and  for  the  gas- 

t 

dynamic  flow  parameters,  i.e.,  for  specific  pg  ,  ,  D^/Dg,  6/D2,  u2, 

the  shape  of  the  leading  edge  of  the  vanes,  etc.  Some  change  in  these 
parameters  (e.g.,  the  thickness  of  the  vanes  6)  compared  with  those 
used  in  the  studies  could  obviously  have  the  result  that  the  influence 
of  z_  would  change  quantitatively  (e.g.,  the  maximum  efficiency  would 
occur  for  some  other  number  of  vanes). 

We  should,  however,  note  that  this  study  was  conducted  for  wheels 
made  according  to  the  generally-accepted  recommendations  for  typical 
wheels  of  stationary  low-output  centrifugal  compressors.  This  fact 
increases  the  value  of  the  obtained  quantitative  relationships. 

In  addition,  the  above  analysis  makes  It  possible  to  draw  certain 
conclusions  as  to  the  influence  of  the  number  of  vanes  on  the  physical 
flow  pattern  In  a  circular  cascade. 

Thus,  analysis  of  the  obtained  experimental  data  makes  it  possible. 
In  many  cases,  to  give  a  qualitative  (and  in  a  number  of  cases  a 
quantitative)  estimate  of  the  influence  of  z_  for  other  design  and 
gasdynamic  parameters  of  the  wheel  and  the  flow. 

In  thus  recommending  that  the  above  data  be  used  In  practice, 
the  authors  at  the  same  time  consider  that  the  above  given  theoretical 
analysis  should  now  be  further  refined,  and  special  experimental 
operations  should  be  undertaken. 

2.  To  investigate  shrouded  wheels  with  normally  backward  bent 
vanes  (in  the  given  case,  S  =  49°),  the  curves  t)*n  and  „  as 

go  clO  ITIclX  IUclX 

functions  of  z  have  a  sloping  nature  for  a  broad  range  of  vane  numbers. 
Therefore  the  maximum  efficiency  of  the  head  coefficient  is  found  for 
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a  broad  range  of  vane  numbers  z  «•  10-20  (for  the  head,  z  *  14-20; 
for  the  efficiency,  z  *<  10-18).  We  should  note  that  the  restriction 
of  the  channel  at  the  inlet  to  the  Investigated  wheels  was  considerable 
(when  z  «  10-20  *  0.78-0. 57)  because  of  the  large  relative  thick¬ 

ness  of  the  milled  vanes  (fi^/D^  =  0.0545)*  These  relations  with 
respect  to  and  6^/D^  are  characteristic  of  low-output  centrifugal 
compressors.  For  wheels  witi.  a  smaller  relative  vane  thickness 
(for  stamped  vanes  61(l  +  (2A/b^))/D^,  where  A  is  the  length  of  the 
flange)  we  can  expect  to  obtain  maximum  efficiency  and  head  coeffi¬ 
cients  for  somewhat  more  vanes. 

Such  an  assumption  can  be  made  with  respect  to  wheels  having  an 

angle  B  which  is  larger  than  that  of  the  investigated  wheels. 
gl 

3.  The  number  of  vanes  greatly  influences  the  value  of  the  zone 
of  stable  and  economic  wheel  operation  regimes  and,  consequently, 
the  width  of  the  acting  zone  of  the  characteristic.  With  an  increase 
in  z,  the  position  of  the  efficiency  maximum  shifts  toward  the  region 
of  lower  flow  rates,  while  the  maximum  head  coefficient  shifts  toward 
the  region  of  lower  flow  rates,  while  the  maximum  head  coefficient 
shifts  toward  the  region  of  higher  flow  rates  and  thus  the  acting 
zone  of  the  characteristic  is  greatly  reduced.  In  this  regard,  to 
increase  the  acting  zone  of  the  characteristic  as  small  a  number  of 
vanes  z_  as  possible  must  be  selected. 

With  large  z,  when  there  is  great  restriction  of  the  flow  by  the 
vanes,  the  acting  zone  of  the  characteristic  can  evidently  be  expanded 
if  we  undercut  every  other  vane,  since  this  should  increase  and, 
consequently,  decrease  . 

•;  4.  The  most  satisfactory  flow  structure  behind  the  wheel,  as 
the  results  of  our  tests  on  the  series  of  wheels  showed,  is  obtained 
for  z  =  10-20. 
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For  wheels  with  z  »  10-20  the  absolute  nonuniformity,  with 
respect  to  the  angle,  of  the  flow  behind  the  wheel  remained  practically 
constant.  In  this  case  the  deviation  of  the  flow  angle  from  the 
average  value  was  +(3-5)°(Aa  -  amax  -  a^n  =  6-10°)  at  a  distance  of 
6.5  ran  behind  the  wheel.  This  should  be  taken  Into  consideration 
when  designing  vaned  diffusers. 

5.  The  angles  of  attack  in  regimes  of  maximum  efficiency  were 
within  the  limits  i±  =  3-7°  when  t±  =  O.78-O.57  (z  =  10-20). 

6.  Comparison  of  the  magnitudes  and  nature  of  the  change  in  the 
and  M"  numbers  (with  a  change  in  tp  ,  z,  t.  )  resulted  In  the  con- 

c  W1  Po  1 

elusion  that  the  number  M  (more  exactly  the  number  M  )  is  not  a 

c  c2 

characteristic  criterion  of  the  operation  of  the  wheel,  but  can  serve 
primarily  as  a  characteristic  of  the  operating  conditions  of  the 
diffuser.  The  criterion  which  determines  the  operation  of  the  wheel 
is  the  number  or  other  numbers  which  are  determined  by  the  inlet 
velocity.  Therefore  there  must  be  theoretical  and  experimental 
studies  to  establish  the  critical  values  of  and  for  character¬ 
istic  designs,  and  these  must  be  used  as  the  criteria  when  selecting 
the  maximum  permissible  peripheral  velocities. 

In  testing  the  indicated  series  of  wheels  for  u^  =  215  m/sec 
in  the  region  z  =  10-18,  where  the  maximum  values  of  max  were 
obtained,  the  value  of  was  within  the  limits  0. 5-0. 53  (Mc  =  0.41- 
0.44). 

Thus,  for  wheels  similar  to  those  investigated  z  =  14-18  should 
be  considered  the  best  number  of  vanes,  from  the  standpoint  of  obtain¬ 
ing  maximum  efficiency  and  head  and  satisfactory  flow  structure  behind 
the  wheel.  If  there  are  rigid  requirements  fbr  the  value  of  the  acting 
zone  of  the  wheel  characteristic,  z  -  10-12  should  be  considered  to  be 
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the  best  number  of  vanes.  For  the  investigated  wheels  the  experimental 


data  show  that  the  formulas  proposed  by  a  number  of  authors  [6,  7,  12, 
13 ]  give  quite  accurate  results  for  finding  the  optimum  number  of 
vanes  with  respect  to  efficiency  and  head.  However,  for  a  proper 
selection  of  the  number  of  vanes,  taking  into  account  the  efficiency, 
the  head,  the  width  of  the  acting  zone  of  the  characteristic,  and  the 
flow  structure  behind  the  wheel,  it  is  necessary  to  use  generalized 
experimental  data. 


Fig.  1.  Geometric  parameters 
of  the  investigated  wheels. 


Fig.  2.  Comparative  character¬ 
istics  of  the  wheels. 


Optimum  Number  of  Vanes  for  the  Investigated  Wheels 
According  to  Various  Recommendations 


Author 

Ponaula 

Opt  laa 

nafcar  af 
van*  a 

Sturt* 

1.  Uttar* 

Tln-§7 

(7,8-6.l).ln  • 

18—14 

m 

0% 

1M  157 
•j  ■  0,35  —  0,45 

V.  r.  Ria 

WKBEStM 

15-16 

161 

■Ul 

C.  Pflaldarar 

The  author  recounendi  that  this  formula 
be  used  for  centrifugal  pu*>e«  In  this 
eese»  k  ■  6.5,  For  centrifugal  eoq?ree* 
sore#  k  ■  11 

22-23 

B 

s.  I.  Taltkla 

wharav,  <  10—  12* 

10-8 

iioj 

'  Ju  I.  Sta panov 

•  *  ~y 

16 

mi 

L.  R.  Voxlta 

*  —  10  +  30Dfc 

wharaO.  1*  In  Batar* 

Raoocaaandad  far  tingla  ahroudad  wheal* 
•f  eantrifugal  eea>r*aaora  of  tranaport 
lnatollatlona 

18 

U«J 

V.  It.  Bdalav 

.sA2slYv* 

* °'07S  din  '  "•  •  , 

19 

\ 

• 

[13J 

' 

•  ^ir‘  ^^«'S(4r) ' . 

Here 

*--h>sr*+K 

S  -  aunfao*  of  van*  waahad  by  goal 

1  -  lanffth  of  davalopatant  tf  vanaj 
■  -  avarag*  valua  of  aagnltudaa  batwaan 
aaatlona  1-1  and  2-2  (a**,  Pig.  2, 
p.  IS). 

40' 


? 


Fig.  J.  Individual  character¬ 
istics  of  wheels  with  various  z’s. 


Fig.  4.  Comparative  character¬ 
istics  of  wheels  for  ri*> _ . 

'ao  max 
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STUDY  OP  THE  EFFICIENCY  OF  CENTRIFUGAL  COMPRESSOR 
WHEELS  WITH  A  TWO-STAGE  VANE  CASCADE 

S.  A.  Anisimov,  et  al. 

At  the  present  time  the  stationary  centrifugal  compressor  produc¬ 
tion  Industry  uses  wheels  with  a  simple  (zg  =  2z^)  and  complex  (zg  > 

>  2z1)  two-stage  vane  cascade  (.1,  2,  3,  4,  5].  The  guiding  principles 
in  this  case  are  basically  that  such  cascades  can  assure  the  required 
flow  guidance  to  the  outlet  from  the  wheel  with  moderate  restriction 
at  the  inlet. 

The  question  of  the  efficiency  and  expediency  of  using  wheels 
with  z2  =  2z^  is  of  great  practical  significance,  since  to  the  present 
time  in  the  literature  there  have  been  no  studies  made  on  the  influ¬ 
ence  of  the  number  of  vanes  and  the  relative  diameter  D^  =  D^/Dg 

(Fig.  l)  on  the  efficiency  of  "compressor-type"  wheels  (B  =  40-50°) 

g2 

with  a  simple  two-stage  cascade. 

These  concepts  also  served  as  the  basis  for  conducting  the  pro¬ 
posed  work,  a  theoretical  and  experimental  study  of  the  influence  of 

z  and  Dk  on  the  operational  efficiency  of  centrifugal  wheels  (0  =  49°) 

”  s2 
with  a  simple  two-stage  vane  cascade.  Certain  data  on  this  problem 

have  been  published  previously  [6]. 
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The  geometric  parameters  of  the  Investigated  wheels  are  given 


in  Fig.  l. 


1.  Theoretical  Study 


As  is  known,  flow  in  the  impeller  is  three-dimensional.  However, 
study  of  the  spatial  flow  in  the  vane  channels  is  extremely  difficult. 
In  the  general  case,  analysis  of  three-dimensional  flow  can  be  replaced 
by  separate  study  of  two  flows  in  the  meridional  and  radial  planes 
with  subsequent  superposition  of  the  flow  patterns  one  upon  the  other. 
To  simplify  the  solution  of  the  problem  it  is  expedient  to  restrict 
ourselves  to  a  study  of  flow  only  in  the  radial  plane,  taking  into 
account  the  change  in  height  of  the  channel  b  in  the  meridional  plane. 
Such  an  assumption  is  possible,  since  usually  the  axial  velocity  com¬ 
ponent  within  the  wheel  is  Insignificant  and  has  no  substantial  effect 
on  the  velocity  profile  In  the  vane  channels. 

On  the  basis  of  these  concepts  we  conducted  an  approximate  theo¬ 
retical  analysis  of  the  Influence  of  z  and  D,.  on  the  operational 
efficiency  of  a  centrifugal  wheel,  with  the  aim  of  obtaining  the 
maximum  head  and  minimum  local  M  numbers  in  the  vane  channels  for  an 
optimum  regime,  and  we  also  estimated  the  infl\ience  of  these  same 
geometric  parameters  on  the  width  of  the  acting  zone  of  the  wheel 
characteristic  and  Its  hydraulic  efficiency. 

The  useful  head  of  the  wheel  can  be  represented  as  follows: 

hn  =  hj*  Afty  —  Lhr, 


where  h^  -  l/g(l-(cr  /u2)cotPg  )u|  is  the  theoretical  head  for  an 

infinite  number  of  vanes;  -  h^  is  the  decrease  in  the  theoretical 
head  in  association  with  a  finite  number  of  vanes  due  to  the  Influence 
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of  the  relative  eddy,  the  curvature  of  the  vanes,  and  the  change  in 
the  flow  structure  in  the  channels;  Ahr  ■  h^  -  h^  is  the  loss  of 
head  due  to  friction  and  vortlcity. 

4 

Let  us  examine  the  dimensionless  value 


The  smaller  4>( z )  the  greater  the  useful  head,  since  the  value  of 
h„,  does  not  depend  on  the  number  of  vanes  or  on  D  .  Consequently, 

•loo  1C 

expressing  the  value  <l>(z)  as  a  function  of  z  and  Dr  we  cam  analyze  the 
influence  of  z  and  on  the  useful  head. 


The  value  f^z )  can  be  represented  as  follows: 


t' ,S  —Af  ,  hr  , 

=  - l-~r 


(2) 


where  p.  -  hp/h^  is  the  coefficient  of  a  decrease  in  the  transmitted 
energy. 

The  value  Ahr  can  be  expressed  as  follows: 


A  hr=Rx 


(3) 


where  Rx  is  the  resistance  arising  with  friction  in  the  vane  channels; 
w^  is  the  mean  relative  velocity  taking  into  account  restriction 
between  sections  1-1  and  2-2  of  the  impeller;  0  is  the  weight  flow 
of  the  gas. 

On  the  basis  of  Kearton's  experimental  data  [7],  let  vis  express 
the  coefficient  for  wheels  with  a  single-stage  cascade  according 
to  Ekkert»s  formula  [4]: 


»*  = 


1 


1  + 


1.5+ I.l^ 


(4) 
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The  expediency  of  using  Formula  (4)  has  been  confirmed  by  our 
experimental  data. 

Determination  of  the  function  fg(z)  is  a  more  complex  problem 
due  to  the  lack  of  sufficient  experimental  data  for  finding  the  resist¬ 
ance  R  .  To  approximately  determine  the  resistance  occurring  with 
turbulent  gas  flow  in  the  vane  channels  of  an  annular  cascade,  let 
us  use,  by  analogy,  the  familiar  relationship  from  the  theory  of 
fluid  flow  in  channels  of  arbitrary  cross  sections 


where  X  is  the  resistance  coefficient;  pm  is  the  mean  density  of  the 
gas  in  the  channel; 


C  . .  o  i ..  i.  \  i  _  fttm  —  ^  "f  bm2  O,  —  O,  _ 

i  «.  U‘„  t  om)  i - ^ 

am  =  ttD^/z  sin  p^  -  6  is  the  width  of  the  vane  channel;  6  is  the 
vane  thickness;  and  1  ■  (D2  -  D^/2  sin  p^  is  the  arbitrary  length 
of  the  channel. 

The  values  of  the  parameters  with  the  subscript  m  were  defined 
as  the  arithmetic  mean  between  sections  1-1  and  2-2. 

The  resistance  coefficient  X,  according  to  recommendations  by 
Balje  [8],  can  be  calculated  from  the  Blasius  formula 


where 


J_  U.3IU 


tt^gm 


Re. 

'  "«i  h^m  * 


Vm 


(6) 


vm  is  the  coefficient  of  kinematic  viscosity.’ 

For  wheels  with  a  single-stage  cascade,  using  Eqs.  (l)-(6)  we 


get,  after  transformations, 


1 


15  -i-i 


i  +■ 


Yw  (»i  +  at)* 


wc«5 


fir) 


2t 


0.810  (Oi~Oi) 


I  _  9,,’oot^,  (  nDm  tin 


(■-*) 


nDm  tin  {hm  \« 
i  —  i*/ 


tin 


ngw  tin  (Wm  -  6*  +  *m» 


n>,  HrJ»»  tin  femfr*, 

Vm  '  tin  fig*,  —  6*  +  bmz 


To  determine  ^(z)  for  wheels  with  a  two-stage  cascade  we  can 
use  Eq.  (7),  leaving  f^z)  unchanged*  and  representing  friction  ar  . 
vorticity  losses  of'  the  entire  wheel  as  the  sum  of  the  losses  in  the 
two  stages 

f2(z)  =  +  f2^z2^‘ 

Figure  2  shows  the  graphic  solution  of  Eq.  (7)  for  the  investi¬ 
gated  type  of  wheels  with  single-  and  two-stage  cascades.  For  wheels 
with  a  single-stage  cascade  the  optimum  number  of  vanes,  with  respect 
to  useful  head,  is  within  the  range  zQpt  h  =  16-25.  Analysis  of  the 
curves  shows  that  with  Increasing  z,  the  optimum  wheel  with  a  two- 
stage  cascade  is  obtained  for  high  D  .  From  the  standpoint  of  obtain- 
ing  minimum  losses  in  useful  head  it  is  expedient  to  use  wheels  with 
a  two-stage  cascade  beginning  with  z  >  18-20. 

We  should  mention  the  restricted  sphere  of  application  of  Eq.  (7) 
because  of  the  approximateness  of  Eq.  (6)  which  defines  the  value  of 
the  resistance  coefficient  X  without  taking  into  account  the  presence 
and  development  of  breakaway  zones  in  the  rotating  vane  channels  of 


*  The  experimental  data  obtained  verify  the  expediency  of  using 
Ekkert*s  equation  [43  to  calculate  the  influence  of  a  finite  number 
of  vanes  of  wheels  with  a  two-stage  cascade,  assuming  z  ■  z2. 
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the  wheel,  and  also  without  taking  into  account  the  effect  of  the 
rotation  and  curvilinearity  of  the  channels.  For  wheels  with  z2  *=  z^ 
the  expediency  of  using  (7 )  is  limited  to  the  region  In  which  local 
diffusivity  and  the  exposure  angles  of  the  channels  are  within  toler¬ 
able  limits  and,  consequently,  the  presence  of  developed  breakaway 
zones  is  not  expected,  i.e.,  a  region  with  a  sufficiently  large  number 
of  vanes.  Below  we  give  an  analysis  which  makes  it  possible  to  estab¬ 
lish  :z  for  which  the  level  of  local  diffusivity  will  be  minimal.  For 
wheels  with  z2  =  2z1  the  restrictions  on  the  equation  for  determining 

X  are  made  more  severe  by  the  fact  that  with  very  high  D  very  short 

/c 

vanes  disrupt  normal  gas  flow  in  the  channels  of  the  second  stand. 

In  addition,  in  the  region  of  very  large  D  it  is  hardly  possible  to 
reliably  determine  the  coefficient  m-  from  known  theoretical  and  empiri¬ 
cal  relationships. 

It  can  be  assumed  that  the  nature  of  the  dependences  «*  f(<p  ) 

-  ^0 
and  ■  f(®  )  will  be  substantially  influenced  by  the  breakaway 

p0 

zones  which  develop  in  the  channels  of  a  centrifugal  wheel;  the  pres¬ 
ence  and  scale  of  these  zones,  in  turn,  depend  mainly  on  the  value 
of  local  diffusivity.  The  value  of  maximum  local  diffusivity  can  be 
estimated  by  the  relation  wloc  max/w2*  Considering  that  the  relative 
velocity  w2  for  a  single  series  of  wheels  with  u2  =  const  changes 
insignificantly,  the  maximum  local  diffusivity  and,  consequently,  the 
degree  of  development  of  breakaway  zones  can  be  estimated  by  the  value 

or  Such  an  approach  to  the  study  of  the  influence 

xoc  max  jloc  max 

of  breakaway  zones  on  the  operational  efficiency  of  a  centrifugal 
wheel  considerably  simplifies  the  theoretical  analysis,  making  it 
possible  to  simultaneously  obtain  a  value  for  the  maximum  local  Mach 
numbers  which  characterize  the  operation  of  the  cascade. 


-49- 


There  exist  various  theoretical  methods  for  calculating  the 
local  relative  velocity  in  the  radial  plane  of  a  wheel  with  a  finite 

number  of  vanes.  However,  most  of  these  methods  are  quite  unwieldy 

/ 

and  when  making  calculations  they  require  use  of  the  method  of  suc¬ 
cessive  approximations.  A  number  of  authors  [1,  2]  have  proposed 
for  this  purpose  a  reasonably  simple  approximate  analytical  dependence 
which  establishes  the  connection  between  the  local  relative  velocity 
and  the  coordinates  of  the  examined  point.  Considering  that  maximum 
local  relative  velocity  is  attained  on  the  nonacting  surface  of  the 

vane,  the  proposed  equation  for  wheels  with  p  <  90°  is  written  as 

g2 

follows : 

%,=Kv{1  ~  2^*)+““. 


where  w^y  is  the  average  relative  velocity  in  the  investigated  sec¬ 
tion,  taking  into  account  restriction  and  bend  of  the  flow  to  aiv  angle 
PgJ  £  ls  the  width  of  the  channel  in  the  radial  plane;  ^  is  the  radius 
of  curvature  of  the  vane;  and  co  is  the  angular  speed  of  rotation  of 
the  wheel. 

Equation  (.8)  is  applicable  for  determination  of  the  local  rela¬ 
tive  velocity  only  in  regimes  of  maximum  airfoil  cascade  flow.  Other¬ 
wise,  due  to  a  change  in  shape  of  the  flow  lines,  there  is  a  consider¬ 
able  change  in  their  equidistance  from  the  center  line  of  the  vane. 

Bearing  in  mind  that 


.  c 

W  =  - a- 

»v  yF 


_G_ 

Y 


and  using  the  expression  for  the  Mach  number  and  Eq.  (8),  we  get, 
finally, 

°n  — i/'T'T  f  K  Yo  r,  'Hi)  ,  «  1  ..  a  nP  sink— ft*  )  (9) 
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A  graphic  solution  of  Eq.  (9)  for  the  investigated  wheels  is 
shown  In  Fig.  3.  An  analysis  of  the  curves  in  this  figure  shows  that 
for  wheels  with  zg  «=  z^  the  optimum  number  of  vanes,  from  the  stand¬ 
point  of  the  minimum  Mach  level  or  minimum  local  diffusivity,  is  within 
the  range  zop^  M  =  11-13*  For  wheels  with  Zg  ■  2z^,  with  increasing 
as.  the  minimum  level  of  the  development  of  local  diffusivity  at  the 
second  stage  is  observed  for  larger  D  .  For  any  ss  the  maximum  value 

of  occurs  when  T5„  =  0.  527 >  l*e.,  at  the  inlet  to  the  first 

np  cz0  K 

stage  of  the  cascade. 

In  two-stage  cascades  of  the  investigated  type,  to  obtain  minimum 
M  numbers  in  the  first  stage  we  must  select  z1  =  10-14,  while  for 
minimum  M  numbers  in  the  second  stage  we  must  assume  D  =  0.7-0. 8. 

K 

Taking  into  account  the  assumptions  used  to  derive  Eq.  (9),  it 

is  to  be  expected  that  the  absolute  values  of  M"  /M  are  approxi- 

wnp  cz0 

mate.  However,  use  of  Eq.  (9)  makes  it  possible  to  quite  simply 
estimate  the  region  of  the  origin  of  the  Intense  development  of  break¬ 
away  zones  In  vane  channels  as  a  function  of  the  geometric  relation¬ 
ships  of  the  wheel  and  the  flow  parameters. 

To  determine  the  Influence  of  the  number  and  length  of  the  vanes 
on  the  hydraulic  efficiency  of  a  wheel  we  can  use  the  functions  f^Cz) 
and  fg(z)  obtained  previously: 

(10) 


The  graphic  solution  of  Eq.  (10)  Is  shown  in  Fig.  3.  Analysis 
of  the  curves  shows  that  with  Increasing  B  ,  when  z  «  const,  the 
hydraulic  efficiency  of  the  wheel  increases.  For  the  same  number  of 
vanes  Zg  the  efficiency  of  a  wheel  with  a  single-stage  cascade  is 
less  than  that  of  a  wheel  with  a  two-stage  cascade.  The  absence  of 
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a  maximum  on  the  x\  curves  is  explained  by  the  above-mentioned 

s* 

restrictions  of  Eq.  (6)  which  does  not  take  into  account  the  presence 

of  developed  breakaway  zones  in  vane  channels  for  wheels  with  relative 

/ 

sparse  cascades  with  a  large  channel  exposure  angle.  Using  previously 
obtained  data  on  the  influence  of  _z  and  on  the  value  of  local 
diffusivity  and,  consequently,  the  level  of  development  of  breakaway 
zones,  we  can  approximately  limit  the  sphere  of  applicability  of  Eq. 
(10)  by  the  curve  M-M.  To  the  left  of  line  M-M  Eq.  (10)  should  yield 
a  qualitatively  incorrect  result. 

Using  simple  physical  concepts  it  is  easy  to  show  that  for 
increased  width  of  the  acting  zone  of  tte  wheel  characteristic  it  is 
necessary  to  decrease  restriction  of  the  flow  at  the  inlet  [9].  This 
assumption  has  been  confirmed  by  experimental  data  obtained.  Conse¬ 
quently,  the  curve  M-M  and  Eq.  (9)  determine  the  optimum  number  of 
vanes  from  the  standpoint  of  attaining  the  maximum  width  of  the  acting 
zone  of  the  wheel  characteristic  for  a  reasonably  high  efficiency 
level. 

It  should  be  mentioned  that  the  proposed  equations  (7),  (9),  and 
(lO)  are  sufficiently  general  and  can  be  used  to  determine  the  influ¬ 
ence  on  the  useful  head,  the  level  of  local  diffusivity,  and  the 
hydraulic  efficiency  of  not  only  z_  and  but  also  of  other  geometric 
relationships  of  the  wheel  and  the  flow  parameters. 

Thus,  on  the  basis  of  this  theoretical  analysis  we  can  assume 
that  wheels  with  two-stage  cascades  can  have  advantages  over  wheels 
with  single-stage  cascades  with  respect  to  efficiency,  head,  and 
width  of  the  acting  zone  of  the  characteristic. 
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2.  Experimental  Study 


The  experimental  part  of  the  work  was  done  on  the  ETsK-1  test 
stand  at  the  Laboratory  of  Compressor  Machines  of  M.  I.  Kalinin  Lenin¬ 
grad  polytechnic  tfnlyerslty  [10 3.  The  diagram  of  the  blading  of  the 
stand  and  the  control  sections  for  measuring  the  flow  parameters  Is 
shown  on  p.  16.  The  wheels  were  tested  for  constant  peripheral 
velocity  Ug  *  214.5  m/sec.  In  all,  25  shrouded  wheels  with  fully 

milled  vanes  were  tested  (P  ~  49°).  Of  these,  8  had  single-stage 

s2 

and  17  had  two-stage  cascades.  The  basic  geometric  parameters  of  the 
wheels  are  shown  in  Pig.  1  and  given  in  Table  1.  In  Table  1,  the 
crosses  designate  wheels  tested  on  the  ETsK-1  stand,  while  the  terms 
"starting"  and  "derived"  cascades  indicate,  respectively,  single-stage 
cascades  with  z  =  z2  and  z  -  z^. 

The  experimental  data  were  measured  and  processed  according  to 
the  Leningrad  Polytechnic  Institute  method  [6,  9,  11 ]  using  small 
pneumometric  devices  developed  at  LPI  [12,  15,  14,  15 J,  and  also  other 
measurement  and  recording  apparatus  [6,  10-163. 

As  criteria  for  estimating  the  efficiency  of  the  wheels  we  used 
parameters  examined  in  detail  in  previous  papers  [9,  11 3.  In  addition 
to  the  parameters  derived  in  these  works,  in  the  present  investigation 
we  used  the  following: 


where 


VWf/ 


(11) 


W. 


sin  flp,  sin  ’ 


(12) 


t  is  the  coefficient  of  flow  restriction  at  the  outlet  from  the  first 

rv 

stage  of  the  cascade  or  at  the  inlet  to  the  second  stage,  determined 
from  z  -  or  z  "  z2^tkII^’ 
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The  values  of  T  and  y  were  determined  approximately,  on  the 
basis  of  the  assumption  of  a  linear  law  of  change  in  temperature  and 
density  of  the  gas  between  sections  1-1  and  3-3,  and  the  presence  of 
the  equalities  TQ  ”  Tl  **  Ti  an{*  *"  ^1* 

Analysis  of  the  individual  characteristics  of  the  wheels  has 
made  it  possible  to  make  the  following  remarks  (Fig.  4): 

1.  The  general  nature  of  the  individual  characteristics  of 
wheels  with  z2  »  2z1  does  not  differ  from  the  corresponding  character¬ 
istics  for  z2  »  z1  [9,  11]. 

2.  The  nature  of  the  curves  ri*-v  =  f(<p  )  and  i/s*  «  f(<p  )  is 

ao  Po  po 

determined  by  the  number  of  vanes  in  the  starting  single-stage  cascade 
z  »  zp  and  the  relative  diameter  13  .  With  an  increase  in  D  for  the 
given  series  of  wheels,  the  curve  becomes  more  inclined.  A  similar 
phenomenon  is  noted  with  decreasing  z  for  D  =  const. 

™  fc 

Such  a  result  is  completely  regular.  An  increase  in  decreases 
the  total  working  surface  of  the  vane 3  and  restriction  at  the  inlet 
and  leads  (when  D  =  l)  to  a  decrease  in  the  total  number  of  vanes  by 

ft 

a  factor  of  two.  Therefore  an  increase  in  and  a  decrease  in  z  of 
the  starting  cascade  should  have  a  qualitative  influence  on  the 
characteristic  of  the  wheel.  This  inter-dependence  of  and  z_  is 
evident  for  other  wheel  characteristics  as  well. 

3.  For  all  the  studied  wheels  T)*^  max  is  within  the  region  i^  = 

«  0-10°,  while  f*.  is  within  the  region  i.  =  15-24°. 

IucLX  J- 

W  K 

increase  with  Increasing  flow  rate.  In  this  case  the  relation  between 
the  values  and  depends  on  the  value  As  a  rule, 

however,  with  large  z  and  very  small  undercuts  we  have  the  inequality 

# 

Kk  >  \l* 


4.  The  number  ^  and  particularly  the  numbers  and  M" 


Analysis  of  the  comparative  characteristics  of  wheels  with  two- 
stage  cascades  makes  it  possible  to  make  the  following  remarks  (Pig.  5)* 

1.  With  increasing  and  decreasing  z  the  zone  of  economic  and 
stable  wheel  operation  increases. 

‘  2.  The  value  of  t}*^  depends  on  T>k  and  z.  The  dependence  max  = 
«  fCS^)  has  an  extremal  value.  For  all  series  of  wheels,  with  increas¬ 
ing  1)  there  is  first  an  increase  in  the  efficiency  to  its  maximum 
value,  then  it  decreases,  and  finally,  when  >  0.8-0. 9,  it  again 
increases  somewhat.  The  nature  of  these  curves  can  be  explained  as 
follows . 

Transition  to  a  two-stage  cascade  due  to  undercutting  at  the 
inlet  of  the  vanes  results  in  decreased  restriction  of  the  flow  (an 
increase  in  t)  and  an  increase  in  the  channel  exposure  angle  v  of  the 
first  stage  compared  with  the  same  sections  of  the  starting  cascade. 

For  the  first  stage  of  vanes,  with  a  further  increase  in  the  values 
of  and  v.  remain  constant,  t  increases,  while  v  decreases  slightly. 
For  rather  large  z_,  there  will  be,  at  the  starting  cascade,  great 
restriction  and  a  small  channel  exposure  angle  at  the  inlet.  A  small 
T1  st  (larse  )  will  have  a  negative  influence  on  the  operation 

of  the  cascade.  The  influence  of  small  st  will,  however,  be  insig¬ 
nificant.  By  undercutting  the  vane  of  such  a  starting  cascade  we 
first  decrease  the  negative  influence  of  great  flow  restriction.  The 
influence  of  undercutting  the  vanes,  however,  on  the  increase  of  the 
angle  v ^  is  less  clearly  expressed.  Therefore,  when  undercutting  the 
vanes  in  the  starting  cascade  with  high  v  we  can  likely  expect  increased 
efficiency.  A  further  increase  in  T5  will  Involve  a  decrease  in  flow 
restriction  and  flow  surface  friction  against  the  vanes.  This  leads 
simultaneously  to  an  increased  load  on  the  vane  apparatus  and. 
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consequently,  to  Increased  local  Mach  numbers  and  local  dlffusivity 
on  the  airfoil,  and  also  to  an  increase  in  the  part  of  the  channel 
having  a  large  exposure  angle.  Because  of  these  factors,  beginning 
at  a  certain  D  there  is  a  reduction  in  efficiency. 

The  negative  influence  of  these  factors  for  large  undercuts,  for 
the  given  wheels  with  >  0.8-0. 9,  is  intensified  by  the  fact  that 
the  short  vanes  of  the  second  stage  of  the  cascade  cease  to  operate 
as  aerodynamic  airfoils  and  are  essentially  additional  resistance. 
Therefore,  with  further  undercutting  of  the  vanes,  right  up  to  com¬ 
plete  shearing  of  the  vanes,  there  is  a  certain  incree.se  in  efficiency. 

These  concepts  make  it  possible  to  assume  that  if  the  starting 

cascade  has  small  z_,  i.e.,  large  st  (small  M^  )  and  large  angle 

1  st 

V1  st'  undercutting  of  the  vanes  may  not  increase  the  efficiency  but, 
conversely,  lower  it.  In  this  case  a  positive  influence  of  the  two- 
stage  cascade  on  wheel  operation  can  be  attained  by  installing  addi¬ 
tional  short  vanes  at  the  outlet.  This  decreases  the  local  dlffusivity 
of  the  channel  at  the  exit  and  evidently  decreases  the  Intensity  of 
the  turbulent  region  on  the  nonacting  surface  of  the  vanes.  This 
analysis  shows  that  the  larger  the  z  of  the  starting  cascade,  the 
larger  will  be  the  values  of  Dr  for  which  the  optimum  operating  regimes 
of  a  two-stage  cascade  will  be  attained.  Experimental  data  have  con¬ 
firmed  this  conclusion. 

3>.  The  head  coefficient  ls  a  function  of  and  £.  With 

a  relatively  small  number  of  vanes  (z  «=  16-24)  for  the  starting  cas¬ 
cade  an  increase  in  Dr  causes  a  decrease  in  With  Increasing 

z  the  region  of  decreasing  decreases. 

The  decrease  In  with  small  z_,  particularly  for  slight  under¬ 
cuts,  Is  evidently  associated  with  a  sharp  (two- fold)  Increase  in  the 
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angle  in  the  cascade  of  the  first  stage  (compared  with  the  starting 
cascade),  while  the  other  factors  are  weakly  expressed.*  We  should 
note  that  with  relatively  small  z_  the  conditions  of  flow  inlet  to  the 
cascade  influence  the  direction  of  the  outlet  velocity  and,  consequently, 
the  head. 

For  large  z,,  when  the  angles  and  small  and  their  influence 
is  not  substantial,  undercutting  of  the  vanes,  which  lowers  the 
restriction  of  the  flow  and  decreases  energy  lost  to  friction  of  the 
gas  against  the  vanes,  can  cause  an  increase  in  the  load  on  the  vane 
apparatus  and  the  diffusivity  of  the  channels  of  the  first  stage.  In 
connection  with  the  above,  the  nature  of  the  change  in  f°r  under¬ 

cutting  of  the  vanes  in  a  two-stage  cascade  will  depend  on  the  rela¬ 
tion  between  these  coefficients  at  the  starting  and  derived  cascades. 


4.  The  values  of  <PpQOpt  and  q^ 


depend  on  15^  as  well  as  on  z. 


With  increasing  D  for  all  wheel  series  there  is  an  increase  in 

/c 

9  and  a  decrease  In  „T,.  The  greatest  change  In  cp^ _ and 


P0°Pt 


rP0cr^_  —  . . - rP0°Pt 

sone  D  SI  O.55-O.62.  For  large  I)  the 

rC  fv 

remain  practically  constant  and  only  when  I)  > 


q>  is  observed  in  the  zone  D  =  0.55-0-62.  For  large  5  the  values 
pQcr  K  K. 

°f  fp0opt  Vr 


>  0.85-0.  9  and  for  small  z_  is  there  again  a  certain  decrease  in  q> 


P0cr 


for  large  D  ,  while  for  large  £  (series  24/12)  there  -is  no  decrease 


whatsoever. 


The  shift  in  the  maxima  of 


''ad 


and  f*  vs.  T)  and  z  can  be  explained 

/C 


by  the  influence  of  the  change  in  the  value  of  and  on  the  ratio 
between  Mj  and  M”  ,  and  and  cp. 

1  -L  1C  K 

The  abrupt  shift  of  the  maximum  efficiency  toward  higher  flow 
rates  when  iL  £  0.58-62  can  evidently  be  explained  by  the  sharply 

rC 


*  We  should  mention  that  the  question  of  the  influence  of  the 
angles  v  on  operation  of  the  cascade  is  a  complex  problem,  one  that 
;*.v.  been  little  studied.  Special  studies  are  required  to  solve  it. 
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decreased  restriction  of  flow  at  the  inlet  and,  consequently,  by  the 


sharp  decrease  in  compared 


with 


*v- 


An  increase  in  in  final  analysis,  decreases  flow  restriction 

/ 

in  a  two-stage  cascade  and  therefore  it  in  principle  should  lower  the 

value  of  cr.  The  greatest  decrease  in  q>p  cr  for  £  0.58-0.62 

is  evidently  explained  in  the  same  way  as  the  increase  in  <pp  Qpt  in 

this  region.  The  decrease  in  ©  for  D  >  O.85-O.9  can  be  associated 

PQCr  K 

with  the  fact  that  further  undercutting  of  the  vanes,  eliminating  the 
inefficient  airfoils  of  the  second  stage  of  the  cascade,  results  in 
improved  flow  aerodynamics  and,  consequently,  an  increase  in  the 
ir'»  <  m"  numbers  in  the  cascade  channels.  Therefore  retention  of 


\  5  \cr 


the  inequality  M",  £  H*'  _  is  now  possible  for  larger  i,  (smaller  <p  ). 

W1  wi  Pq 

With  increasing  z  and,  consequently,  decreasing  load  on  the  vanes  and 

angle  v,  the' influence  of  greatly  shortened  vanes  on  the  flow  aero¬ 
dynamics  decreases.  Therefore  there  is  less  of  a  decrease  in  <pp  cr 

in  the  region  D  >  0. 85-O. 9,  while  for  the  series  of  wheels  beginning 
/c 

with  24/12  there  is  none  whatsoever. 

Analysis  of  the  curves  k  =  f(D  )  shows  that  the  width  of  the 

y  K 

acting  zone  of  the  characteristic  depends  on  D  and  on  z.  In  the 

IQ  "" 

region  D  £  0. 58-O.62  there  is  a  sharp  expansion  of  the  acting  zone 
of  the  wheel  characteristic.  A  further  increase  in  Dk  does  not  result 
in  a  substantial  change  in  k„,  and  only  when  5..  >  O.85-O.9  for  small 

j  K 

£  is  there  again  observed  an  increase  in  ky.  There  is  less  expansion 
of  the  acting  zone  of  the  characteristic  in  this  range  of  changes  of 
I}^  with  increasing  z,  while  for  the  series  of  wheels  beginning  with 
24/12  there  is  no  expansion.  The  nature  of  the  change  in  ky  Vs.  15^ 
and  z  is  determined  by  the  corresponding  shift  in  the  maximum  effi* 
ciency  and  head  coefficient. 
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5-  The  relative  flow  nonunifomity  with  respect  to  energy  (the 

coefficients  k£  and  k^  i)  does  not  depend  on  and  z. 

The  relative  flow  nonuniformity  with  respect  to  angle  and  flow 

rate  (the  coefficients  and  kr)  depends  on  I)  .  With  an  increase  in 

to  O.76-O.8  the  coefficients  kQ  and  k^  remain  practically  constant. 

A  further  increase  in  D  leads  to  an  increase  in  k  and  k„.  This 

tc  a  g 

latter  can  evidently  be  explained  by  the  above-mentioned  negative 

influence  of  the  greatly  shortened  vanes  of  the  second  stage  of  the 

cascade.  Removal  of  these  vane3  again  decreases  k  and  k_.  As  is  to 

be  expected,  the  surge  on  the  curves  of  k  and  kr  for  D  >  O.76-O.8 

ct  y  /c 

corresponds  to  a  certain  decrease  efficiency  in  this  zone. 

The  shape  of  the  curve  of  the  absolute  flow  nonuniformity  with 
respect  to  the  angle  Aa  =  f(D  )  is  basically  identical  with  that  of 

/C 

the  curve  k_  =  f(B  ),  and  can  be  explained  by  the  sake  reasons. 

CZ  K 

6.  The  numbers  M  and  M,f  depend  on  z  and  also  5  for  q*-*  . 

—  fC  ctO  lUaJC 

The  numbers  and  ,  with  increasing  increase  somewhat.  This 

is  due  to  the  fact  that  an  Increase  in  results  in  an  increase  in 
9?  opt  and,  consequently,  in  w^  and  w^  and,  respectively,  Mw  and  . 
With  increasing  z  cp  decreases,  which  lowers  M  .  The  value 

-  TP0°Pt 

of  with  increased  is  greatly  influence  by  the  increased  flow 
restriction  which  increases  w^,  and  therefore  increases. 

J>.  Comparison  of  the  Results  of  a  Theoretical  Analysis 
with  Experimental  Investigations 


Wheels  with  a  Single-Stage  Cascade.  The  proposed  equation  (7) 
for  estimating  the  Influence  of  various  geometric  criteria  of  a  wheel 
and; the  flow  parameters  on  the  value  of  the  useful  head  can  be  used 
also  to  determine  the  optimum  number  of  vanes.  The  graphic  solution 


of  this  equation,  given  in  Pig.  2,  shows  that  the  optimum  number  of 
vanes  from  the  point  of  view  of  attaining  maximum  useful  heads  for 
the  investigated  wheels  is  in  the  region  zopt  ^  -  16-25 .  From  experi¬ 
mental  studies  zQpt  ^  ->  14-20.  Thus,  the  experimental  data  obtained 
eonfirm  the  expediency  of  using  theoretical  dependence  (7). 

The  possibility  of  using  Eq.  (9)  for  determining  the  optimum 
number  of  vanes  from  the  standpoint  of  achieving  the  maximum  width 
of  the  acting  zone  of  the  wheel  characteristic  with  a  sufficiently 
high  efficiency  has  been  confirmed  by  experimental  data  obtained. 

For  the  investigated  wheels  z  ^  M  =  10-14  vanes.  According  to  experi¬ 
mental  results  zQpt  M  =  11-15*  In  this  case,  compared  with  a  wheel 
with  z2  -  z1  =  16  ky  increases  2-  to  5-fold,  while  the  efficiency  only 
decreases  1. 5-2.0^. 

Wheels  with  a  Two-Stage  Cascade.  The  proposed  equation  (7),  from 
a  qualitative  standpoint,  correctly  evaluates  the  influence  of  z  and 
on  the  losses  of  useful  head  in  a  wheel  with  a  two-stage  cascade. 
Analysis  of  this  equation  shows  that  with  increasing  D  the  optimum 
number  of  vanes  zQpt  h  increases.  From  the  point  of  view  of  achieving 
minimum  losses  of  useful  head,  use  of  the  two-stage  cascade  is  Justi¬ 
fied  only  when  z2  £  18-20.  This  agrees  with  experimental  data  obtained, 
in  accordance  with  which  z2opt  h  >  24. 

The  use  of  Eq.  (7)  to  calculate  zQpt  h  in  the  region  of  large 
leads  to  quantitatively  incorrect  results  because  of  the  mentioned 
restrictions  on  dependence  (6)  due  to  the  negative  influence  of  the 
greatly  shortened  vanes  of  the  second  stage. 

Comparison  of  theoretical  and  experimental  data  makes  it  possible 
to  confirm  that  Eq.  (9)  can  be  used  to  determine  the  optimum  values 
of  z  and  for  wheels  with  a  two-stage  cascade  from  the  point  of  view 
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of  attaining  maximum  width  of  the  acting  zone  of  the  characteristic 
with  reasonably  high  efficiency.  The  theoretical  and  experimental 
curves  In  the  region  D  <  0.8,  whosn  in  Fig.  6,  practically  coincide. 

Shrouded  wheels  having  a  single-stage  cascade  of  normally  back¬ 
ward  bent  vanes  with  z2  "  z±  =  16-24  are  widely  used  in  stationary 
centrifugal  compressors.  The  conversion  from  single-stage  cascades 
of  such  wheels  to  two».stage  cascades  with  z2  =  2z1  makes  it  possible 
to  increase  the  efficiency  and  expand  the  acting  zone  of  the  wheel 
characteristic  (Fig.  6).  True,  with  relatively  small  £  (in  our  case 
z2  =  16  and  20),  this  Involves  a  certain  reduction  in  t)*^  and  ky> 
as  well  as  with  respect  to  f*. 

The  relative  increase  (positive  or  negative)  in  the  values  of 
max'  ky'  and  ^opt  for  two-stage  cascades  of  each  series, 

compared  with  their  starting  single-stage  cascades  and  with  an  optimum 
single-stage  cascade  having  z  =  16,*  is  shown  in  Table  2. 

Our  experimental  study  has  made  it  possible  to  make  a  qualitative 
and  quantitative  estimate  of  the  influence  of  £  and  on  the  operating 
efficiency  of  a  wheel  with  a  two-stage  cascade.  The  quantitative  data 
obtained,  strictly  speaking,  are  valid  only  for  specific  wheel  and 
flow  parameters,  i.e.,  for  specific  pgl,  p  ,  pg2,  D^/D,,,  6/D2,  shape 
of  the  leading  edge  of  the  vanes,  etc.  A  change  in  these  parameters 
compared  with  those  used  in  the  experiments  might  have  the  result  that 
the  influence  of  z  and  D  would  change  quantitatively,  to  a  certain 
extent.  It  should  be  noted,  however,  that  our  study  was  made  for 
wheels  constructed  according  to  generally  accepted  recommendations 


•  *  Investigations  of  the  LPI  [8]  showed  that  the  optimum  single- 
stage  cascade,  from  the  standpoint  of  maximum  efficiency  and  head  for 
the  given  wheel  geometry,  is  obtained  when  z2  =  =  16. 
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for  typical  wheels  of  stationary  low-output  centrifugal  compressors. 

Use  of  the  obtained  theoretical  dependences  makes  possible  a 
qualitative  and  quantitative  analysis  of  cascades  for  various  geometric 
dimensions  of  a  wheel  and  distinct  flow  parameters. 


Fig.  1.  Diagram  of  the  investigated 
wheels. 


Fig.  2.  Determination 
of  zopt  h  of  wheels  with 

single-  and  two-stage 
cascades. 
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Fig.  ;5.  The  effect  of  z  and  DR 

on  the  level  of  local  M  numbers, 
losses  In  useful  head,  and 
hydraulic  efficiency. 


Fig.  5»  Comparative  characteristics  of 
wheels  of  the  series  20/10  ( ^amax  **  22°). 


Fig.  6.  Comparative  characteristic  of  wheels 
with  starting  and  derived  single-stage,  and 
optimum  two-stage,  cascades. 


TABLE  2 

Comparison  of  Wheels  with  Single-  and  Two-Stage  Cascades 


«./«« 

®«opt 

Jonparison  with  stantinj 

».*«.*!•  i 

. 

% 

^ ad 

A** 

'ad 

16,8 

0,59 

+4 

-3.2 

-H 

-3.2 

20/10 

0,65 

+6 

-3.9 

+4 

—2.4 

.24/12 

0.78 

+9 

—4 

+3.5 

— 4.7 

28/14 

0J9 

+  11, 

+6.5 

+4 

+2.3 

32/16 

,0.8 

— 

— 

+  1.0 

+4 

i 
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CALCULATION  OP  VARIABLE  REGIMES  OF  CASCADE  PLOW 


V.  A.  Zysin 

When  investigating  one-dimensional  flow  and,  in  particular,  when 
analyzing  variable  regimes  of  operation  of  a  turbine  stage,  great 
importance  may  be  attached  to  the  clearness  and  relative  simplicity 
of  the  calculation  methods  used.  This  is  particularly  important  in 
the  case  of  supercritical  flow  regimes. 

Here  it  might  be  of  great  help  to  introduce  an  i-S  diagram  with 
a  grid  of  the  dynamic  characteristics  —  lines  of  constant  specific 
flow  rates  (Fanno  lines)  [4,  5)*  This  would  make  it  possible,  e.g. , 
to  directly  determine  the  state  of  the  supersonic  flow  behind  the 
turbine  impeller  if  there  is  overexpansion  relative  to  the  given 
counterpressure,  to  detect  the  possibility  of  flow  choking,  and  to 
solve  a  number  of  calculation  problems. 

However,  for  practical  application,  it  is  necessary  that  the 
diagram  be  generalized.  Here  it  is  desirable  that  the  coordinates  of 
the  dimensionless  parameters  to  be  determined  and  the  selected  scale 
be  most  adaptable  to  the  usual  methods  of  turbine  design. 

'•  In  the  present  article  we  propose  using  a  generalized  diagram, 
similar  in  structure  to  the  ordinary  i-S  diagram,  for  a  calculation 
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analysis  of  the  variable  regimes  of  turbine  cascade  operation. 

Calculations  using  this  method  aim  to  explain  the  most  general 
regularities.  Therefore,  all  processes  are  schematized  to  the  utmost: 
we  use  the  stream  theory,  assume  that  the  velocity  and  flow  factors 
are  identical,  do  not  take  into  account  ejection,  leakages,  etc. 

In  light  of  the  above,  the  methods  herein  proposed  are  not  intended 
to  replace  profound  gasdynamic  calculations  which  strive  to  take  into 
account  to  the  fullest  possible  extent  the  influence  of  various  factors. 

1.  The  Generalized  Dynamic  Diagram 


Figure  1  shows  a  generalized  dynamic  diagram  constructed  for  an 
ideal  gas  with  an  isentropic  exponent  k  =  l.J. 

The  diagram  itself  and  the  method  for  constructing  it  have  been 
described  previously  [lj.  Therefore,  let  us  limit  ourselves  to  a 
brief  examination  of  it. 

The  diagram  was  constructed  for  isoenergetic  flow  with  irrevers¬ 
ible  losses.  Consequently, 

=  i  +  const,  1  (i) 

ds  >  0, 

where  1q  is  the  stagnation  enthalpy,  A  =  1/427  is  the  thermal  equivalent 
of  work,  £  is  the  velocity,  and  £  is  the  entropy. 

Along  the  abscissa  is  plotted  dimensionless  entropy: 


2 


At 

Tr 


In— t*— 
Wi-o 


In  Pt, 


(2) 


where  pQ  is  the  stagnation  pressure  for  the  given  entropy,  (Pq)z  „  0 
is  the  stagnation  pressure  for  z  =  0,  £  is  the  dimensional  static 
pressure,  P  -  p/(pQ)z  Q  is  dimensionless  pressure,  and  R  is  the  gas 
constant. 

For  convenience  in  using  the  diagram  the  abscissa  contains  a 
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uniform  scale  of  dimensionless  pressures. 

Along  the  ordinate  (from  top  to  bottom)  is  plotted  the  dimension¬ 
less  velocity 


x 

fl  CVTt 


(3) 


where  C  =  ^\/2g  k/(k  +  1)  R  is  a  constant,  a*  is  critical  velocity, 
and  Tq  is  the  absolute  stagnation  temperature. 

The  diagram  contains  a  grid  of  dimensionless  isobars  P  =  const 
and  a  grid  of  constant  reduced  flow  rates  (dimensionless  Fanno  lines) 

Here  y  is  the  specific  weight,  and  (cy)*  _  Q  is  the  critical  specific 
flow  rate  (mass  flow  rate)  at  zero  entropy. 

From  the  upper  left-hand  corner  of  the  diagram  are  drawn  flow- 
lines  with  constant  velocity- loss  coefficients  9  or  f. 

If,  parallel  to  the  ordinate,  we  plot  the  vertical  scales  of  M, 
T/Tq,  and  the  other  gasdynamic  functions,  then  when  using  the  diagram 
it  is  not  necessary  to  resort  to  gasdynamic  tables. 

In  addition,  on  the  diagram  we  could  plot  grids  of  the  dynamic 
characteristics  which  would  make  it  possible  to  calculate  the  com¬ 
pression  shocks.  In  Fig.  1,  for  simplicity,  these  additional  con¬ 
structions  are  omitted. 


2.  A  Cascade  for  Subsonic  Velocities 

Figure  2  shows,  in  X-z  coordinates,  the  process  in  a  turbine 
;ascade  for  subcritical  flow  regimes. 

When  using  one  of  the  usual  methods  for  calculating  a  turbine 
for;  a  variable  regime,  the  use  of  a  diagram  can  facilitate  determina¬ 
tion  of  the  state  of  the  flow  behind  the  cascade.  If  the  parameters 
of  the  incoming  flow  are  given. 
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Actually,  knowing  the  dimensional  relative  velocity  of  the  incom¬ 


ing  flow  £  and  the  critical  velocity  a*,  we  can  determine  the  corre¬ 


sponding  dimensionless  velocity  X  -  c/a*.  Prom  this,  on  the  X-z 

c  c 


diagram,  we  find  the  point  a  which  corresponds  to  the  state  of  the 
flow  entering  the  casoade. 

Let  us  calculate  the  reduced  flow  rate  on  leaving  the  cascade: 


<D  _  <p  8lna  h 
•  4  lina,  17* 


(5) 


where  h  is  the  height  of  the  vanes  of  the  preceding  cascade  or  the 
height  of  the  guide  channel,  h.  is  the  height  of  the  vanes  at  the 

C 

exit  from  the  calculated  cascade,  a  is  the  angle  formed  by  the  cascade 
axis  and  the  vector  of  the  incoming  velocity,  is  the  angle  formed 
by  the  cascade  axis  and  the  exit  axis  of  the  vane  channel;  ®a  is  the 
reduced  flow  rate  of  the  incoming  flow,  and  <!>e  is  the  reduced  flow 
rate  on  emerging  from  the  cascade. 

Having  determined  from  certain  concepts  the  velocity  coefficient 
for  the  cascade  q>,  let  us  find,  at  the  intersection  of  the  lines  <I>  = 

«=  const  and  <p,  rf/  =  const,  the  point  b  which  corresponds  to  the  flow 
state  at  the  exit  section.  From  this  we  determine  directly  the  dimen¬ 
sionless  exit  velocity  X  ,  the  corresponding  dimensionless  pressure 

C1 

P  ,  and  the  ideal  exit  velocity  X  . 
e  cQ 

When  using  the  generalized  diagram  it  is  easy  to  detect  certain 
phenomena  which  can  occur  in  variable  operation  regimes. 

For  example,  it  may  be  that  the  line  <l>e  corresponding  to  the 
reduced  flow  rate  at  the  impeller  exit  (Fig.  j)  does  not  intersect 
the  line  <p,  V'*  This  indicates  choking,  during  which  the  exit  section 
of  the  cascade  cannot  handle  the  given  flow.  From  the  diagram  we  can 
immediately  determine  the  maximum  flow  which  the  impeller  can  handle 
for  the  given  parameters  of  total  stagnation  ahead  of  it.  This  flow 
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o.s  deterained  by  the  reduced  flow  4^,  the  line  of  which  touches  the 
curve  i>,  at  point  m. 

3.  A  Cascade  for  Supersonic  Flow  Velocities 

'  Let  us  examine  a  cascade  consisting  of  a  Laval  nozzle.  When 
analyzing  the  process  of  expansion  in  nozzles  we  should  take  into 
account  the  difference  in  losses  in  the  convergent  and  divergent  parts 
of  the  nozzle.  Often  we  can  consider  that  all  losses  are  concentrated 
in  the  divergent  part  [3].  Then  point  k.  on  the  X-z  diagram  (Fig.  4) 
corresponds  to  the  state  of  the  flow  in  the  sonic  throat.  Assuming 
r,hat  the  flow  fills  the  entire  exit  section  of  the  nozzle,  let  us 
;:etermine  the  reduced  flow  in  this  section 

«e  -  Wfe' 

where  f|n1r|  is  the  are  of  the  sonic  throat,  and  f  is  the  area  of  the 
exit  section. 

Similarly,  let  us  find  the  point  b  which  corresponds  to  the  stat^ 

of  the  flow  in  the  exit  section  at  the  intersection  of  the  lines  q>  = 

«  const  and  4>  -  const, 
e 

If  losses  in  the  convergent  part  cannot  be  neglected,  i.e.,  the 
flow  factor  x  <  the  reduced  flow  in  the  sonic  throat 

*  X* 
m 

It  is  appropriate  to  note  that  in  this  case  the  sonic  throat  is 
characterized  on  the  diagram  by  the  point  m  which  has  the  corresponding 
dimensionless  velocity  ^  -  1.* 

.■■■■■■  <i— ■ 

*  From  the  diagram  we  can  immediately  see  that  according  to  the 
stream  theory  the  velocity  is  subsonic  in  the  sonic  throat  of  an 
actual  Laval  nozzle. 
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The  throat  is  defined  fcy  the  point  k'  on  the  intersection  of  the 
line  <p  -  const  (e  for  the  convergent  part  of  the  nozzle)  with  the 
line  X  =  1. 

The  reduced  flow  in  the  exit  section 

^e  "  ^mln^e** 

The  position  of  the  throat  (in  the  divergent  part  of  the  nozzle) 
is  determined  from  the  relationship 

Vfmln  ‘ 

where  <t»K  is  the  reduced  flow  at  the  point  k’ . 

Further,  there  are  possible  two  basic  regimes,  depending  on  one 
of  the  following  conditions: 

1)  the  relative  rate  of  flow  into  the  next  cascade  is  supersonic, 
or 

2)  the  relative  rate  of  flow  into  the  next  cascade  is  subsonic. 

In  the  first  case  we  can  assume  that  the  static  pressure  in  the 

axial  clearance  does  not  change,  and  the  next  cascade  is  calculated 
in  the  usual  manner. 

If  as  a  result  of  this  calculation  it  appears  that  the  state  of 
the  flow  behind  the  cascade  is  characterized  by  a  point  on  the  diagram 
in  the  region  X  <  1,  we  should  be  on  the  lookout  for  choking  in  the 
given  stage. 

In  the  second  case,  when  the  rate  of  flow  into  the  next  stage, 
obtained  from  the  velocity  triangle,  is  subsonic,  pressure  in  the 
axial  clearance  in  the  general  case  will  differ  from  the  pressure  in 
the  nozzle  exit  section,  and  will  depend  on  the  amount  of  pressure 
drop  in  subsequent  stages,  and  on  the  counterpressure. 
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If  the  pressure  in  the  axial  clearance  pn  is  defined  in  some  way 
or  other,  the  cascade  operating  in  the  region  of  supersonic  velocities 
is  calculated  as  a  function  of  the  sign  of  the  inequality 

^  Pn> 

where  pe  is  the  static  pressure  in  the  exit  section  of  the  cascade. 

If  p  >  p„  *  p  ,  there  should  he  divergence  in  the  transverse 
e  n  ru. 

section.  Let  point  b  (Fig.  4)  define  the  state  of  the  flow  in  the 
nozzle  exit  section.  This  section  has  the  corresponding  dimension¬ 
less  pressure  p  and  the  dimensionless  velocity  X.  (in  the  case  of 
e  c± 

convergent  nozzles  point  b  can  lie  only  on  the  line  X  =  l). 

Let  us  determine  the  dimensionless  pressure  in  the  axial  clearance 


pn1  *  pn/p0< 

where  pQ  is  the  total  pressure  (static  and  dynamic)  ahead  of  the 
nozzles. 

The  state  of  the  flow  behind  the  transverse  section  is  defined 

by  the  point  n^  at  the  intersection  of  the  isobar  Pn  with  the  line 

<p,  ^  =  const.  The  dimensionless  velocity  X  _  >  X.  and  the  reduced 

clnl  C1 

flow  4>  correspond  to  the  point  n, .  From  the  ratio  <X>  /<D  we  can 

-L  e 

find  the  divergence  in  the  transverse  section  by  one  of  the  familiar 
methods.  If  we  consider  that  there  are  no  losses  in  the  transverse 
section,  point  b1  on  the  vertical  dropped  from  point  b  onto  the  isobar 
Pr  corresponds  to  the  state  of  the  flow  behind  the  section. 

When  the  pressure  in  the  cascade  exit  section  is  lower  than  that 

in  the  axial  clearance,  strong  shock  waves  are  unavoidable. 

# 

,  If  we  assume  that  the  flow  Is  not  separated  from  the  walls  at 

; 

the  end  of  the  vane  channel,  the  state  behind  the  exit  section  is 

defined  by  the  intersection  of  the  line  4>  with  the  dimensionless 

© 
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isobar  corresponding  to  the  given  pressure  in  the  axial  gap.  Figure 
4  shows  cases  when,  depending  on  the  counterpressure,  the  terminal 
velocity  remains  supersonic  (the  point  rig)  or  becomes  subsonic  (the 
point  n^).  Actually,  within  the  nozzle  and  in  the  exit  stream,  there 
ocours  a  system  of  compression  shocks  accompanied  by  separation  of 
flow  from  the  walls.  Naturally,  it  is  impossible  to  analyze  such 
phenomena  on  the  basis  of  a  model  of  one-dimensional  flow  used  in  the 
given  calculation  method.  At  present  it  is  most  correct  to  use  exper¬ 
imental  data. 

However,  we  can  state  that  the  calculation  method  given  here 
gives  maximum  losses.  Actually,  if  there  is  flor  separation  from  the 
walls  (or  vice  versa)  on  intersection  of  the  shock  front  in  the  free 
stream,  the  current  density  increases.  Consequently,  the  new  reduced 
flow  <I>  is  greater  than  the  reduced  flow  9  which  corresponds  to  flow 

©  v 

without  separation.  But  now  expansion  to  the  given  dimensionless 
pressure  p^  or  p^  is  given  not  by  the  points  ng  and  n^,  but  by  the 
points  n£  and  n^  respectively.  Figure  4  shows  directly  that  in  this 
case  the  value  of  the  terminal  entropy  decreases,  while  the  velocity 
increases. 


4 .  Calculation  of  a  Cr  cade  for  ChOiij-ng 

Above  w^  examined  the  case  when  the  cascade  cannot  pass  tr e  give 
flow  due  to  insufficient  area  of  the  exit  sections. 

It  has  been  shown  experimentally  that  often  flow  choking  is 
caused  by  the  angle  of  attack  at  the  cascade  inlet  [2].  Actually, 
we  know  that  a  decrease  in  current  density  corresponds  to  a  flow 
shock  into  the  concave  part  of  the  vane,  while  an  increase  in  curreir. 
density  corresponds  to  a  flow  shock  into  the  convex  part.  In  this 
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latter  case  there  Is  the  danger  of  flow  choking  by  the  leading  edges 
of  the  cascade. 

If  the  relative  velocity  of  the  Incoming  undisturbed  flow  has  the 
corresponding  reduced  flow  <t>00,  the  reduced  rate  of  flow  Into  the  cas¬ 
cade 

®ln  =  sinCPi  -  6)/sin  (6) 

where  6  Is  the  angle  of  attack,  considered  positive  for  a  shock  Into 
the  concave  surface;  and  is  the  inlet  angle  of  the  vane  apparatus. 

Obviously,  even  in  the  ideal  case  it  is  necessary  that 

*in  <  1‘ 

From  this  condition  we  can  determine  the  negative  angle  of  attac  c 
for  which  there  will  be  no  flow  choking  on  entry  into  the  cascade. 
Considering  the  finite  thickness  of  the  leading  edges  and  the  addi¬ 
tional  irreversible  losses,  we  should  consider  that  for  subsonic 
velocity  of  the  incoming  flow  we  must  assume  the  probability  of  chok¬ 
ing  at  <I>e  >  O.975,  to  which  correspond  (at  k  =  l.jj)  the  values  of  the 
dimensionless  velocity  Xin  >  O.85.  In  the  case  of  slight  overlappings., 
for  thick  leading  edges,  choking  occurs  even  earlier. 

When  the  incoming  flow  is  supersonic,  the  danger  of  choking  exist 
continually  with  negative  angles  of  attack.* 

It  may  be  that  the  Bection  limiting  the  flow  is  located  in  the 
vane  channel,  if  in  the  latter  the  cross  section  is  minimum. 

The  cascade  in  this  case  can  be  preliminarily  checked  as  follows. 

We  construct  a  graph  on  which  along  the  abscissa  is  plotted  the 

/  *  For  example,  experiments  conducted  at  Leningrad  Polytechnic 
Institute  and  the  Central  Scientific  Research  Institute  for  Boilers 
and  Turbines  revealed  choking  in  all  cases  when  the  supersonic  incom¬ 
ing  flow  had  negative  angles  of  attack  [2]. 
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chord  of  the  airfoil  b,  and  along  the  ordinate  is  plotted  the  reduced 
flow,  using  as  unity  the  critical  flow  on  entry  into  the  cascade 
(PiS-  5). 

,  t 

The  critical  reduced  flow  at  the  cascade  exit 

%  -  V 

where:  P„  is  the  dimensionless  pressure  at  point  c  (Pig.  2). 

If  we  assume  that  there  is  uniform  loss  in  head  along  the  chord 
of  airfoil  b,  dz/db  4=  const.  The  change  in  critical  reduced  flow  in 
the  vane  channel  is  thus  depicted  by  the  straight  line  A-B  (Pig.  5)« 

On  this  same  graph  we  can  plot  the  curve  of  the  change  in  reduced 
flow  in  the  channel: 

0  =4.  =  f  /f  , 
n  a  a'  n* 

where  <J>  is  the  reduced  flow  in  the  given  section,  and  f0  Is  the  area 
n  0  a 

of  the  corresponding  section,  determined  from  the  drawing. 

If  the  line  <&n  intersects  line  A-B  which  corresponds  to  the 
critical  flows  (segment  C'-D’  in  Pig.  5)>  we  should  suspect  that  the 
vane  channel  will  not  be  able  to  handle  the  given  flow. 
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Fig.  1.  Generalized  diagram 
for  velocity  (x  -  c/a+)  minus 
entropy  (z  -  As/AR).  Abscissa  — 
scale  of  dimensionless  pressure. 
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Fig.  4.  The  process  in 
Laval  nozzles. 


Fig.  2.  Representation  of 
the  process  in  an  airfoil 
cascade  on  a  X-z  dir.gram 
for  X  <  1. 
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Fig.  5.  Case  in  which  the 
exit  section  will  not  handle 
the  given  flow. 


Fig.  5-  Checking  the  vane 
channel  for  choking. 


FTD-TT-62-1816-1+2+4 


-79- 


REFERENCES 


1.  V.  A.  Zysin.  Grafoanaliticheskly  metod  rascheta  odnomernykh 
potokov  szhimayemoy  zhidkos ti.  Trudy  LIP,  No.  2,  1954. 

/ 

2.  V.  A.  Zyaln,  et  al.  Teploenergetika,  No.  2,  1955* 

3.  I.  I.  Kirillov,  and  S.  A.  Kant  or.  Teoriya  i  konstruktaiya 
parovykh  turbin,  Mashgiz,  1947. 

4.  A.  Stodola.  Darapf  und  Gaaturbinan,  VI  Auflage,  1924. 

5*  J.  Kostin.  J.  M.  E.  Proc.,  Vol.  159*  1948. 


FTD-TT-6l-l8l6-l+2+4 


■80- 


DISTRIBUTION  LIST 


DEPARTMENT  OF  DEFENSE  Nr.  Copies 


HEADQUARTERS  USAF 


AFCIN-3D2  1 

ARL  (ARB)  1 


OTHER  AGENCIES 

CIA  ) 

NSA  6 

DIA  9 

AID  2 

OTS  2 

AEC  2 

PWS  1 

NASA  1 

ARM!  3 

NAVI  3 

RAND  1 

POE  12 

NAFBC  1 


MAJOR  AIR  COMMANDS  Nr.  Cipio 


AT  SC 

SCFDD  1 
ASTIA  25 
TDBTL  5 
TDBDP  5 

AEDC  (AZT)  1 
SSD  (SSF)  2 
BSD  (BS?)  1 
AFFTC  (FTT)  1 
AFSWC  (SWF)  1 


FTD-TT' 


62-1816/1+244 


81 


